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Abstract 

We review in detail the Hamiltonian dynamics for constrained systems. Emphasis is put on 
the total Hamiltonian system rather than on the extended Hamiltonian system. We provide 
a systematic analysis of (global and local) symmetries in total Hamiltonian systems. In par- 
ticular, in analogue to total Hamiltonians, we introduce the notion of total Noether charges. 
Grassmannian degrees of freedom are also addressed in details. 
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1 Introduction 



Symmetries have always been a determinant guide for the understanding of Nature, because they 
seem to enable the simultaneous concretization of the two ideals underlying the scientific quest: 
simplicity and beauty.^ The symmetry principles have been essential for the development of 
modern physics, e.g. in the birth of both relativity theories or in the building of the standard 
model. The importance of symmetries has been recognized since the very beginning of scientific 
inquiry but mankind waited until the twentieth century for a new paradigm to emerge: the gauge 
symmetry principle.^ The aphorism "symmetry dictates interaction" can be considered as the cor- 
nerstone of modem theoretical physics. Both in classical general relativity and in quantum field 
theory, the gauge symmetries are the deep geometrical foundations of fundamental interactions. 
Indeed, gauge symmetries determine the terms which may appear in the action. Nevertheless, 
some qualifications need to be made because, unfortunately, the symmetries rarely fix uniquely 
the interactions although this dream underlies most unification models. Even though, at first 
sight, gauge transformations could have been naively dismissed as auxiliary -if not irrelevant- 
tools since they are in some sense "unphysical," they actually proved to be almost unavoidable! 
For instance, from a field theoretical point of view the light-cone formulation is perfectly con- 
sistent by itself, but it is extremely convenient to introduce spurious unphysical (in other words, 
"gauge") degrees of freedom in order to write down Lagrangians for massless particles which 
are manifestly local and covariant under Lorentz transformations. Another example is general 
relativity where the decisive role played by the requirement of covariance under the diffeomor- 
phisms does not need to be stressed, even though a superficial glance at this issue would dismiss 
this requirement as irrelevant since any theory can be formulated independently of the coordinate 
system by introducing an affine connection.^ 

Like for every deep and fundamental concept in physics, gauge symmetries exhibit many faces 
and can be approached in different ways. The investigations of Dirac on the Hamiltonian formu- 
lation of gravity opened a new door for entering into the world of gauge theories. As explained in 
his seminal works [3], the presence of gauge symmetries in the Lagrangian framework implies, 
from a Hamiltonian point of view, the existence of "constraints" on the phase space variables. 
Conversely, the study of constraints in the Hamiltonian framework may serve as a path leading 
towards some understanding of gauge symmetries. The present lecture notes are intended to be 

^One may recommend the collection of inspiring lectures given by Chandrasekhar or the celebrated book of Weyl 

on this topic [ 1 ] . 

^The many developments of this crucial chapter in the history of physics are very well summarized in the book [2]. 

'a student, scared by some of the conceptual subtleties arising from the gauge symmetry principle, could find 
some recomfort in the following surprising anecdote: during his quest for a reconciliation between gravity and 
relativity, Einstein himself initially argued that the equations of motion for the metric must not be diffeomorphism 
covariant! The "physical" bases of this wrong initial requirement were related to the subtle issues mentioned above. 
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a self-contained introduction to the Hamiltonian formulation of systems with constraints. Since 
the seminal investigations of Dirac, the development of this topic has been so dramatic that we 
would not pretend to be complete. At most, we do hope that these notes might be useful to new- 
comers searching for a pedestrian and concrete approach on the interplay between Lagrangian 
vs Hamiltonian systems from the point of view of gauge symmetries vs constraints. The main 
particularity of the present notes is that the (rigid and gauge) symmetries and their associated 
Noether charges are discussed with many details. For instance, the various possible definitions 
according to the choice of formalisms (Lagrangian, total or extended Hamiltonian) are introduced 
and compared with each other. Their explicit relationship is provided, due to its importance for 
applications. Another original feature is that the fermionic case is included in the presentation 
from the very beginning. This case is so relevant in physics that we found better to discuss the 
general case immediately in order to allow a uniform treatment of all physical cases, rather than 
devote later a specific section to this 'particular' case. We also insert all the details and proofs of 
the properties of supermatrices that are used in this text. More generally, pretty much all results 
presented here are given with their proofs, in order to be entirely self-contained. Nevertheless, 
we have not aimed at complete mathematical rigour (in the sense that all symbols written in the 
text are supposed to exist under suitable regularity conditions and that the formal manipulations 
they are subject to, are allowed). The major emphasis is on the classical level, though we provide 
some flavour of the quantization process at the end of these notes, for both first and second class 
constraints. 

In contradistinction to most of the fundamental textbooks on the subject, such as [4-7], we 
focus on the total Hamiltonian instead of the extended Hamiltonian."^ On the one hand, the main 
advantage of this choice is that the dynamics determined by the former is always equivalent to the 
Lagrangian dynamics. On the other hand, its drawback is precisely that the primary constraints 
play a privileged role, while such a distinction is not relevant from a purely Hamiltonian perspec- 
tive. Of course, the evolution of the physical quantities (i.e. the observables) through the dynam- 
ics of either the Lagrangian, the total Hamiltonian or the extended Hamiltonian always agrees. 
Therefore, the preference between total and extended Hamiltonian is somehow 'philosophical', 
in the sense that it reflects the opinion whether, respectively, the Lagrangian formulation is more 
fundamental than the Hamiltonian one, or the contrary. Mathematically, one may argue that none 
of these opinions is more valid than the others because some Lagrangian systems do not allow 
an Hamiltonian formulation, and conversely. (Both types of examples are reviewed in this text.) 
Physically, the quantization process seems to rely heavily on the Hamiltonian formulation, even 
if Feynmann's path integral could plead in favour of the Lagrangian as well. Still, a rigorous 

"^Of course, these textbooks do include very detailed discussions on the total Hamiltonian formalism, we only 
mean that it is not their chief emphasis and guideline. Of course, the formalism of constraints is discussed in many 
other textbooks, e.g. [8,9]. 
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definition of the path integral measure, etc, seems to be more natural in terms of the phase space. 
Although the total Hamiltonian formulation is emphasized here because our perspective is on 
the relationship with the Lagrangian formulation, we will adopt an 'oecumenic' attitude by dis- 
cussing both approaches. 

The plan of these notes is as follows: The Lagrangian dynamics is reviewed in Section 2 
under complete generality, scrutinizing on various aspects of the symmetries of the action prin- 
ciple which are not always addressed in details (higher derivatives, finite versus infinitesimal 
transformation, invertibility, etc). The canonical Hamiltonian formalism for a dynamical system 
with constraints is reviewed in Section 3 and the conditions of the equivalence between the La- 
grangian and the Hamiltonian formalism are mentioned. In the section 4, the total and extended 
Hamiltonian dynamics are introduced together with the distinct types of constraints: primary 
or secondary, first or second class. Although these distinctions become somewhat irrelevant at 
some deeper level from the Hamiltonian point of view, they are very important when addressing 
the quantization process or the concrete relation between the Lagrangian and the Hamiltonian 
formulations. The symmetries of dynamical systems and their associated conserved charges are 
discussed thoroughly in Section 5 from several perspectives. The equivalences between the many 
approaches are not shown through general theorems but through direct computations in order to 
provide for the reader a concrete grasp of the formalism. The Dirac quantization, suitable for 
second class constraints and reviewed in Section 6, has a rather straight interpretation: eliminate 
the spurious degrees of freedom by making use of the Dirac bracket. The main drawback of this 
method is that, in general, one is unable to compute the Dirac bracket explicitly. While the BRST 
quantization, suitable for first class constraints and presented in Section 7, is more subtle concep- 
tually (because it is formulated as a cohomological problem) and technically (because many fields 
such as ghosts, etc, have to be added) it possesses at least one great virtue: if the Hamiltonian con- 
straints (or the Lagrangian gauge symmetries) have been entirely determined, then this method 
can always be settled concretely in order to write down the gauge-fixed path integral, even if 
the BRST cohomology group cannot be computed explicitly. At the end come some appendices: 
some proofs of various properties are placed in the appendix A in order to lighten the core of the 
text. A rigorous treatment of the fermionic variables is provided in Appendix B through a review 
of the Grassmann algebras, while all the necessary definitions and basic properties of superma- 
trices are provided in Appendix C. The proofs of some propositions on the canonical forms of 
supermatrices are presented in details in the appendix D. Finally, the appendix E is devoted to 
a very simple and illustrative example of the general discussion contained in the body of these 
notes. We advise the reader to progressively go through this example, while (s)he goes through 
the general material in the core of the text. 

These notes are an expanded version of some lectures given by JHP at Sogang University 
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during the years 2007 and 2008 . 
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2 Lagrangian dynamics: symmetry and Grassmann variables 
2.1 Euler-Lagrange equations 

We first consider a generic Lagrangian depending on M variables,^ q^, 1 < A < J\f, their time 
derivatives, and time is allowed to appear explicitly, 

J0.{qn,t), (2.1) 

where 

-j n = 0,l,2,3,-.-. (2.2) 

The 5^'s form the coordinates of the so called "jet space". Note that some of the variables can 
be fermionic.^ In our conventions, unless explicitly mentioned, all derivatives act from the left to 
the right, 

dFl^d£^ r_i^#A(#F+#A)lZ (2 3) 

where is the Z2 -grading of the 2jV-dimensional tangent space with coordinates (g^^, g^), 

„ _ J for bosonic A 

~ \ 1 for fermionic A ' ^ 

From the variation of the action 

S[q^]^ J C{qn,t)dt, (2.5) 
and up to the boundary terms, one obtains 

J dt SC{qn, t) = j dt Sq' 



E' '''' 



_n=0 

SO that the corresponding Euler-Lagrange equations 

SC{qm,t) 



dtj dq^ 



C{qm,t), (2.6) 



6q^ 



= 0, (2.7) 



^For some issues in the continuous limit Af ^ 00, see e.g. [9]. 

^One of the only prerequisite of these lectures is that the reader is supposed to be familiar with graded algebras 
and related super objects. A good self-contained introduction to supersynnmetry is [10]. 
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are given by acting on the Lagrangian with a Unear differential operator called the Euler-Lagrange 
operator/ 



n=0 ^ ^ 



n 



In the jet space, the submanifold defined by Eq.(2.7) is called the "stationary (hyper)surface". 
We have introduced the symbol = which will stand, from now on, for 'equal on the stationary 
surface' or, equivalently, 'equal modulo the Lagrangian equations of motion (2.7)'. Remark that 
the Euler-Lagrange operator (2.8) is not a derivation, i.e. it does not obey to the Leibniz rule. 
It is useful to understand that in the jet space the total derivative d/di is defined as 

d d ^ A d , 

n=0 

d d d d d d d d d 

This implies that the Euler-Lagrange equations of a total derivative term are identically vanish- 
ing, 

0. (2.11) 



SO that 



Sq^ V 

The converse is also true, therefore 

j^{qn,t)^0 <^ F{qr,,t)^ . (2.12) 

We present a proof of these statements in Appendix A. Note that a generalization of Eq.(2.12) 
holds for field theories {J\f = oo) too, in which case it is referred to as "algebraic Poincare 
lemma". ^ 



2.2 Symmetries in the Lagrangian formalism 

In general, a symmetry of the action involves a certain change of variables, 

/ AQn,t), (2.13) 

'For the fermionic degrees of freedom, there arises a subtle point which we discuss in Appendix B. 
^For more details on this lennma, on jet space, etc, see e.g. the section 4 of [11] and references therein. 
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which may explicitly depend on the q^'s and the time t . It corresponds to a symmetry of the 
action if the Lagrangian is invariant under the transformation up to total derivative terms. 



where 



d 



d 



m=0 



(2.14) 



(2.15) 



Surely this imposes nontrivial conditions on the form of g(g„, t) in terms of the Lagrangian, 



One useful identity for an arbitrary function F is^, for m > 1, 



dqn 



dq 



~B 



n=0 ^ / \ in / n=0 

which gives the following algebraic identity for any change of variables, 

n / oo c\~B 



5C{qut) 



oo 



5q^ 

Infinitesimally, q^ 



n=0 
A 



dt) 



\m=0 



dqn 9qg 



n=0 



dV fdq^Sjr,{qi,t) 



dt) \dq^ 5q^ 



qo + Sqoiqn, t), Eq.(2.17) implies 



5q' 



,5 



(2.16) 



(2.17) 



(2.18) 



smce 



r SC{qi,t) _ 5C{qi,t) \ _ ^ SC{qi,t) _ SC{qi,t) 



V 6q' 

oo 



Sq' 



\ 6q^ 



n=0 



dq^ 5C{quf) 
dt) \ dq^ Sq^ 



5q' 
6C{qi,t) 



Sq 



(2.19) 



Eq.(2.17) will be used later in order to show that the Eler-Lagrange equations are preserved by 
symmetries of the action, a fact which is natural to expect but is non trivial to prove. 



'See Eq.(A.lO) for a proof. 
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Now assuming the symmetry (2.13), acting with the Euler-Lagrange operator on both sides of 
(2.14), using (2.11) and (2.17), we get 

5C{qut) ^ 5C{qut) ^ ^ /_^y (dq^ 5C{qut) \ 

Sq^ Sq^ ^ V dt) \dq^ 5q^ ) ' ^ " 

Before we discuss the generic cases, we first consider the simple case where g(f (g, t) depends 
on q and t only, being independent of qn (n > 1) and invertible i.e. det{dq/dq)j^ 0. We have 

5C{qi,t) _ dq^iq,t) 6C{qi,t) ^^^^^ 
Sq^ dq^ dq^ ' 

using (2.20) with q and q exchanged. Hence, if q{t) is a solution of the equations of motion, then 
so is q{q, t), where q{q, t) depends on q and t only, i.e. not on the derivatives g„ (n > 1). 

Now, for the generic cases where q{qm t) depends on the g^'s (n > 0): if there exists an in- 
verse map q{qn, t) - most likely depending on the infinite set of variables'*^ g„ (n = 0, 1, 2, • • •) - 
then the (inverse relation of) Eq.(2.20) indeed shows that if q{t) is a solution of the equations of 
motion, then so is q{qn, t). 



2.2.1 Invertible transformations 

The existence of an inverse map is always guaranteed when there exists a corresponding infinites- 
imal transformation, 

q^^q^ + Sq^, 5q^ ^ f^{q^,t) . (2.22) 

Consequently, the infinitesimal transformation of the coordinates of the jet space are given by 
Qn + (^?n ' where 

Sq^ = {^Jf^i^rn.t) - + E ^m^B^ f^iqm.t) . (2.23) 

More explicitly, we define an exponential map with a real parameter s, 

f °° d \ 

q^{s, qn, t) = exp Is ^ t)-^ j / , g^(0, g„, t) = q^ . (2.24) 



^"For example, the usual translational symmetry reads 

n>0 ' n>0 
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From (2.10) we first note the commutativity property. 



and hence. 



^ j ?m, ^) = exp \ s ^ j ■ (2-26) 



The main claim is then. 



ds 



i=0 



From (2.26), the first equality in (2.27) is obvious. The derivation of the other relation is carried 
in Eq.(A.12) of the appendix. The equation (2.27) implies that the following differential operator 
is 's' -independent, 

E f^i^rn.t)^ = E f^i^rn.t)-^ • (2.28) 
Z=0 1=0 

Using the above identities, it is straightforward to obtain the explicit inverse map, q"^ — > q"^, 

qn, t) - exp f-s/f (g^, t)^) . (2.29) 



2.2.2 Local symmetries 

In the case of a "local symmetry", namely if the transformation involves some arbitrary time 
dependent functions a'^{t) as 



q^(qn,t,a{t)y (2.30) 



it is possible to have 'different' solutions by varying the arbitrary functions ai{t) , even though 
one starts from the same initial data qn{to) . However, one may consider that the Lagrangian 
alone dictates the whole dynamics of the given system (and nothing else) and, furthermore, that 
the dynamics is deterministic (i.e. there is a unique solution to the Cauchy problem). As long 
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as one takes this viewpoint for granted, then one must regard different trajectories as the same 
physical state. 

Namely any local symmetry must be a 'gauge' (i.e. unphysical) symmetry. In mathematical 
terms, a physical state is given by an equivalence class for which the local symmetry defines the 
equivalence relation.^^ Obviously, the presence of the gauge symmetries complicates the correct 
counting of the number of physical degrees of freedom (especially if the gauge symmetries are 
not independent, etc). An "observable" quantity is a function on the space of physical states, 
hence it must be gauge invariant. We will turn back to these issues in the Hamiltonian context. 



2.3 Second order field equations 

Henceforth, we focus on the standard Lagrangian, L{q,q,t), which depends on q^, q"^, t only, 
and derive some algebraic identities for the later use. 

When the infinitesimal symmetry transformation, q"^ — > q"^ + Sq'^{q,q,t), also depends on 
g^, g^, t only, we have^^ 

Thus, the function 6K{q,q,t) in 

6L = ^5K (2.32) 
dt 



^^A complete and detailed treatment of the gauge invariance of an action can be found in the chapter 3 of [7]. 

'^The Hamiltonian framework enables a precise and 'algorithmic' computation of the number of degrees of free- 
dom, which leads to a precise and completely general criterion for counting the physical degrees of freedom directly 
from the form of gauge transformations in the Lagrangian formaUsm itself [12]. This rigorous treatment clarifies the 
origin of some maxims from the physicist folklore (such as "gauge shoots twice") and thereby provides a supple- 
mentary argument in favor of the fruitful interplay between Hamiltonian and Lagrangian formalisms. 

'^Actually, in the case of a regular second-order Lagrangian, this can be assumed without loss of generality, as 
explained in the exercise 3.8 of the book [7]. 
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must depend on q"^, q"^, t only too, as 

A 7- x^A dLiq,q,t) A dL{q,q,t) 

^ ^^^ dL{q,q,t) ^ (^nd{5q^ ^ -sdi^ ^ d{5q^\ dL{q,q,t) ^^ ^^^ 

dq"^ \ dq^ dq^ dt J dq^ 

dt dq^ dq^ dt 

This implies 



(2.34) 



d{5q^)dL(q.q,t) _ d{6K) 
dq^ dq^ dq^ 
and 

By taking the partial derivative of the latter equation with respect to q^ , making use firstly of 
Eq.(2.31) and secondly of Eq.(2.34), we get 



QqA QqB dq^dq^ dq^dq^ \ dq^ dq'-^dq^ dq^dt J dq^ 



d{5K) ^ ^b 9^{^K) ^ d\6K) ^ ^c9{Sq^) d'^L 



dq^ dq^dq"^ dq^dt dq^ dq^dq^ 

dq"^ \ dq^ 9t J \ dq^ dif ) dif dq^dqf 



(2.36) 

From Eq.(2.34) we find that the coefficient of q'^ must vanish, which implies an integrability 
condition: 



dq^ dq^dq^ dq^ dq^dq^ 



'_^^#.#c A (d{5q^dL\ _ d_ fd{Sq^dL\ ^ 
' ^ dq^\ dqC dq^ ) dq^ \ dq"^ dq^ J 



(2.37) 
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The transformation of the 'momenta' dL/dq^ is given by 



/ dL{q,q,t) 
V dq^ 



-Sq' 



QqBQqA 



QqBQqA 



d(5K) _ d{5q^) dL d{6q^ 
QqA QqA QqB QqA 



(2.38) 



dL 



dt \dq^ J dq 



dL 



where Eqs. (2.36) and (2.37) have been used in the derivation at the second Une. Finally, for the 
transformation of the Hamiltonian we get 



.A9L{q,q,t) J, . ,.\ .^"9(^9^) 

L(g,g,i)) = g 



dq' 



dL 



dt \dq^ J dq 



dL 



+ 



d{5q' 
~~dt 



dL 



dq 



d{5K) 



dt 
(2.39) 

These relations will be used later in Section 5.1 where we analyze the symmetries in the Hamil- 
tonian formalism. 



The corresponding Noether charge is defined by 



Up to the Euler-Lagrange equation. 



dL d / dL 



(2.40) 



dq^ dt \dq^ J 



(2.41) 



the Noether charge is conserved. 



dQ , A 



dt 



d f dL\ dL 



dt \dq^J dq 



0, 



(2.42) 



due to Eq.(2.32). Further discussions on symmetries in the Lagrangian dynamics are carried out 
in Section 5.1. 
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3 From Lagrangian to Hamiltonian and vice versa 



3.1 Canonical momenta 

Given a standard Lagrangian L{q, q, t) , depending on M bosonic or fermionic variables (with 
1 < A < M) and their first time derivatives (and, possibly, on time as as well), the equations of 
motion read 



dq^ dt \dq^ J ^ dq^dq^'^'^ dq^ dq^ dtdq^ ' ^^'^^ 

If the jV X jV supermatrix, d^L/dq^dq"^, is nondegenerate, all the g^'s are uniquely determined 
by q and q. Namely all the variables are completely determined by the initial data, and also all 
the 'velocities' may be expressed in terms of q and the canonical momenta '■— dL/dq^ . 
Henceforth, we focus on the degenerate case, 

In a large class of examples, it may still possible that all the variables are uniquely determined 
from the initial data through the equations of motion, e.g. for a Lagrangian L{q, g, t) = uabQ^Q^— 
V{q, t), linear in the 'g^', and where the constant graded symmetric matrix, uJ[ab} = ^ab, is non- 
degenerate. The (anti)symmetrization has weight one, i.e. 

i^[AB} \ {i^AB - . (3.3) 

We will not attempt to analyze and classify all cases here in the Lagrangian formalism, but we 
will do so in the Hamiltonian formalism later. 

3.2 Primary constraints 

Now, let us start from the expressions for the M momenta in terms of g, g, t, 

PA^^^^^§^^fA{q,q,t) A^l,2,...,Af, (3.4) 

and try to invert the map in order to express the velocities g^ in terms of g, t and the momenta p. 

We first consider a bosonic system having bosonic variables only. If one of the J\f momenta 
(3.4), say pa, depends nontrivially on a certain velocity, say g", then this velocity can be expressed 
in terms of pa and the remaining velocities, collectively denoted by q"^, as well as q"^ and t . Then, 
substituting the expression 

g«-/i«(g^,p„,g-,t), (3.5) 
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into the other momenta than pa , collectively denoted by , we get 

Pm = 5'm(?^ , Pa , g"" , . (3.6) 

This procedure can be repeated until the expressions for the momenta pm do not depend on any 
of the velocities. Performing the procedure step by step a finite number of times, we get finally 



r^h^{q^,Pa,q^,t), 



Pm = fmiq^ ,Pa,t), 



where the sets of momenta and velocities split into two disjoint groups each, 

PA^{Pa,Pm) ■■ {a}U{m} = {1,2,- ■■,J\f} , {a}n{m} 



0, 



(3.7) 



(3.8) 



= g'^) : {a} U {m} = {1, 2, • • • , jV} , {a} n {m} = , 

and in particular our procedure defines a one-to-one correspondence of {a} <-> {a}.^"^ In words, 
on one side the velocities have hatted indices, on the other side the momenta have unhatted 
indices. For the velocities, the rh's correspond to the velocities which remain independent while 
the a's correspond to the velocities which are determined in terms of the former velocities and 
momenta. For the momenta, the situation is opposite: the a's correspond to the ones which are 
independent while the m's correspond to the momenta which are expressed in terms of the latter 
momenta. 

Notice that Pmiq, q) = fm (^q^, Pa{q, q),t^ are identities on the tangent space of coordinates 

(g^, q^) and that there exists an invertible map between the momenta Pa and the velocities 
(keeping g^, and t fixed), 

qa^h%q^,Pa,f\t) ^ Pa = faiq"" :q^:t). (3.9) 

The latter means that 

'dpa{q, q,t)' 



det 



dq^ 



(3.10) 



thus the rank of dpA/dq^ is equal^^ to the dimension of the set {a}. Also, one may say that there 



is a one-to-one map 



16 



Pa 



(3.11) 



'''They are distinguished because pa is not necessarily the conjugate momentum of g". 

' ""This fact is an alternative starting point for getting the constraints and the decomposition of the indices in disjoint 
set. We prefered to provide a concrete explanation of the result (3.7) instead of a slightly more abstract one in terms 
of the rank of the ( super) Jacobian via the implicit function theorem. 

'^See (3.34) for a more general result. 
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Now we return to the generic systems having both bosons and fermions. In contrast to the 
bosonic system, the above procedure which expresses the velocities in terms of the momenta may 
not work even if the momenta depend on velocities nontrivially, mainly due to the non-existence 
of an inverse for any fermionic variable. We consider an example, 

L = leex, (3.12) 

which gives pe = i Ox and Px = i 00- When 6 is fermionic while x is bosonic (which is the case 
with the usual notations) none of the expressions can be inverted. Furthermore, the corresponding 
hamiltonian reads H — i 99x — L which again cannot be reexpressed by the momenta and coor- 
dinates only.^^ If 9 were bosonic, then H — —ipoPx9~^. In the present paper we do not consider 
this case. We always assume that when the the momenta depend on velocities nontrivially, one 
can always obtain the inverse function until we achieve the expression (3.7). Namely we will 
restrict'^ our analysis to systems of bosons and fermions, where the 2J\f coordinates {pa, q^) of 
the phase space (1 < ^, -B < AA), are subject to M. functionnally independent constraints, 

0m(p,g,^) = O, l<m<M. (3.13) 

The M. primary constraints are independent in the sense that the following A4 vectors are linearly 
independent, 

S,^™ =(f=,f=,-,|^) , l<m<M, (3.14) 

V \dpi dp2 dpu J V 

where |y means that the left-hand-side is evaluated on the hypersurface defined by the system 
(3.13), after taking the partial derivatives. The implicit function theorem actually ensures that 
it is possible to solve Eq.(3.13) for M. of the momentas, as in Eq.(3.7). Also, the M. primary 
constraints naturally define a {2^ — A^) -dimensional hypersurface V in the phase space, called 
the "primary constraint (hyper)surface", 

l^ = {(p,g)|0m(p,g,i) = O, l<m<M}. (3.15) 

'^One possible way to circumvent the obstacle is to employ explicitly the Grassmann algebra basis and work 
strictly with the real number coefficients, namely \_q'^\j , [pa]j , as discussed in Appendix B. One can then apply 
the above procedure in the bosonic system without any problem until one gets a similar expression to (3.7). However, 
the corresponding Hamiltonian dynamics for all the coefficients, especially the Poisson bracket, will have to be de- 
composed into a complicated expression, and this will not be done here. We will always assume that the expressions 
of the constraints do not need any explicit use of the Grassmann algebra basis. 

^^If the constraints are not independent from each other, they are said to be "reducible." This more general case is 
treated in the section 1.3.4 of [7]. 
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We furthermore assume that at fixed time t all the constraints can in principle be solved to express 
any point on V by 2A^ — M. independent variables (1 < i < 2J\f — M), 



V={{p,q)=f{x,t)} 



(3.16) 



which provide a local coordinate chart on V (the time dependence is due to the fact that the 
primary constraints may depend explicitly on time). On the other hand, M. constraints of all ^^'s 
can be taken as coordinates for the linearly independent directions to V in the full phase space. 
The entire 2j\/'-dimensional phase space has then two sets of coordinate charts. 



Note that there may exist some freedom in choosing different sets of the independent momenta 
{pa} from the M primary constraints.'^ 

For an arbitrary function g, t) on the 2jV-dimensional phase space, with the coordinate 
system [x, (f)) in (3.17), we define F{x, 0, t) :— F{p, q, t) and a set of M. functions, F'^{p, q, t), 
by 



F{p, q, t) = F{x, cf>, t) = F{x, 0, t) + q, t) = F{p, q, t)\y + q, t) , (3.18) 



where F{p, q, t) \y = F{x, 0, t). In other words, |y means that we substitute (p, q) by its expres- 
sion in terms of (x, t) on the stationary surface V . 

Finally note already that when we study the Hamiltonian dynamics, there can appear more 
constraints, namely the "secondary constraints". In this case, all the constraints will define a 
smaller hypersurface, V cV,m the phase space. 

3.3 Prior to the Hamiltonian formulation: change of variables 

In this subsection, instead of (g^, q^) we regard {q"', q^,Pa, q^) in Eqs.(3.7) as the independent 
variables, and discuss briefly the time evolution of them. The Lagrangian equations of motion are 

^^For more details on the regularity conditions and the properties they imply, some of which are used here, the 
reader is referred to the subsection 1.1.2 of [7]. Notice also that more general regularity conditions (e.g. where 
the momenta do not play a distinguished role, or where the constraints are not assumed to be independent) can be 
defined, as is done in the first chapter of [7]. 



{pa , q^) , 



1<A,B<N 



(3.17) 




l<i<2M - M, l<m<M. 
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equivalent to 



dpa f dL{q,q,t) , 



dt V 

tU 

(3.20) 



9* = /i"(«^,P6,<j™,t) 



dt \ dq" 



provided 



Essentially these equations lead to a set of M. algebraic relations on {q^^Pa, q™') by substituting 
(3.19) and (3.21) into (3.20):20 

dq-^dq^^\dq'^J.,^,,dpadt \dq"^).,^,,' ^ 



where m runs from 1 to . Now these ^A algebraic relations can be thought as constraints for 
the M. variables g™. Such constraints fix some of the g™'s but may leave others as completely 
free parameters. Once all the ^"^'s are determined as functions of other variables or as free param- 
eters, the time evolution of the remaining (here, taken to be independent) variables (pa, q"',q^) 
follows from (3. 19) and (3.21). However, the constraints (3.22) are in general nonlinear in cf^ and 
so they are difficult to solve. Below, we move to the Hamiltonian formalism where the indepen- 
dent variables are {q^iPa, q^) rather than {q'^^Pa, One advantage is that the corresponding 
constraints will be linear in q'^ so that we can do a more explicit analysis (see Subsection 4.3). 



3.4 From Lagrangian to Hamiltonian 

Suppose that a given Lagrangian leads to M primary constraints, say (3.7): 

0m(g^, PB, t) := Pm - /m(g^, Pa: t) . (3.23) 

^"in terms of the Hamiltonian i/(<7'^, Pa, t) and the Pois son bracket [ , }p,b. defined later (in Eq.(3. 24) and (4.1) 
respectively), the equation (3.22) can be reexpressed in a compact form: 

H{q,pa,t) + q"'{p„-f„) , Pm - fmj ^ ^ + =0> 

where fm= fm{q,Pa, t) and the explicit velocities are taken to be constant with respect to the phase space deriva- 
tives of the Poisson bracket. Anticipating a bit, one may realize that, in such a way, there are A4 linear equations 
(3.29) for the M. variables q"^ rather than the M. algebraic nonlinear equations (3.22) for g"*. 
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Replacing by h^{q,Pa, cp, t) in (3.23) we take again {q^,Pa, t) as the independent vari- 
ables. We write the "canonical Hamiltonian" as 

H{q^,Pa,t) -.^q^pA- L{q,q,t) 

^ h%q^,Pa,q'^,t)pa + q'^Prh- L {q^,h%q'',pa,q'^,t), q^) (3.24) 

= q''Pa + q"'fm{<l'^,Pa,t) - L{q^,h''{q^,Pa,q"\t), g™) , 

where we made use of Eq.(3.7). Since there may exist some freedom in choosing different sets of 
the independent H — M. momenta, the Hamiltonian is not uniquely specified, in general, from a 
given Lagrangian, but depends on this choice. However, on V, these Hamiltonians are all equal. 

From the fact that 

dH _ dh^ _ dh^ dL _ dL _ 

one can see that the canonical Hamiltonian is indeed a function of q^, Pa and t only, i.e. it is 
independent of q^ . Further direct calculations can lead to 



9H{q,Pb,t) ^ ^_-^.#^.a _ r_-^Wa#mqmd^ 

dpa dpa ' 

dpn dpn 



(3.26) 



dH{q,pb,t) _ fdL{q,q,t)\^ ^ ,,^90, 



dq^ \ dq^ A. ^ ^ dq^ 



where all the velocities are to be understood as functions of {q^,pa, q"^,t) by the substitution 
= h"'{q,pa, q"', f). The first two equations are easily obtained by making use of the Legendre 
transform philosophy, i.e. the canonical Hamiltonian does not really depend on the velocities. 
Concretely, it is enough to perform the partial differentiation only of the momenta in the term 
q'^PA in order to compute the right-hand-side of the first lines from (3.26). In a unified manner, 
any velocity can be expressed as a function of {q,pa,q^,t) 

q^{q,Pa,r,t) (-!)#- ^^^1^ + . (3.27) 

OPA OPA 
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Now this foraiula suggests that we can take not only (g, p^, g'", t) but, alternatively, (g, Pa,c["^, t) 
as independent variables. As follows from Subsection 3.3, the Hamiltonian dynamics is consistent 
with the Euler-Lagrangian equations only if 



dt 

d/m 
dt 



dH{q,p^,t) 
dq"- 

dH{q,pb,t) 
dq"^ 



+ 



dfn 



(3.28) 



(3.29) 



where we made use of the definition (3.24) of the canonical Hamiltonian and of the constraint 
(3.23). The equation (3.29) leads to Ai linear constraints on the A4 variables g"*. Once we know 
the complete solution of g™, as we will do in Sec. 4. 3, the other equations (3.27) and (3.28) deter- 
mine the time evolution of the remaining variables Pa-,(f- 



For an equivalent but more compact description of the Hamiltonian dynamics for the variables 
(g^, Pa, g"*, t), we introduce the "total Hamiltonian" defined by 



HTiq"", pB,u^,t):= Hiq^, Pa,t) + <Prn{q^, pb, t) v!^ , 



(3.30) 



where Ht is indeed a function on the 'total' phase space with coordinates (g^,PA) but demand 
that 

OHt _ dL{q,q) 



T 



dp^ 



{-l)*-q\ 



V 



dq' 



V 



dq' 



Combining (3.23) and (3.31) we identify 



hence dmU 



m^"^ = q"^4>m ■ We also have 



H' 



T\v 



H\ 



V ' 



(3.31) 



(3.32) 



(3.33) 



V 



In summary we note that there exist two one-to-one maps: 

/ g^ \ / g^ \ 



Pb 



Pb 



(3.34) 
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The Lagrangian dynamics is equivalent to the Hamiltonian one with the primary constraints. In 
the former system, the dynamical variables are by definition {q^, q^}, while in the latter the set 
of independent variables can be chosen to be {q^,Pb, u'^} for convenience. 

As we will shortly show in Subsection 4.3, in the Hamiltonian dynamics there is a systematic 
way of identifying the variables {u^}. Once the most general solution u^{q,p,t) preserving 
the primary constraints is obtained, as in (4.45), the Hamiltonian dynamics determines the time 
evolution of the other variables {q^,Ph} , with the restriction on V. After the generalization to 
the 'total Hamiltonian system', it can govern the dynamics of the 2iV-dimensional whole phase 
space, with variables {q^,PB} , free from any restriction. 



3.5 From Hamiltonian to Lagrangian 

We start^^ with a given Hamiltonian H(p, q) on the 2A/'-dimensional phase space of coordinates 
{PA: q^) , and M independent arbitrary but fixed primary constraints, 0m (p, 5, ^) = , which can 
be solved as pm — fm{(l^,Pb, t) to express M. of the momenta in terms of the positions and the 
N — M. other momenta. There can be some freedom in choosing different sets of the independent 
momenta, {pa}- The relevant phase space reduces to a (2A/' — 7W) -dimensional hyperspace V, as 
in (3.15), which will be further restricted to its sub-manifold V <Z V if there occurs 'secondary 
constraints' (see Subsection 4.3). 

Introducing Ai new variables vJ^ , we define the total Hamiltonian, 

M 



Ht(p, q, u, t) = H{p, q,t) + ^ (t)m{p, q, t)u"' . 

m=l 

The action principle is derived from the action 

S[q^,PB,u"']^ J dt{pAq^-HT). 



(3.35) 



(3.36) 



For instance, this leads to 

q^(q,Pi,u\t) 



dp. 



■1) 



V 



OH ^ (90m ^m 
dpA dpA 



(3.37) 



V 



where \v means that we substitute Pm by p^ — fm{q^-, Ph-, t) after taking the partial derivatives. 



^'This subsection is provided for completeness, and may be skipped at the first reading. 
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We assume that there exists an inverse map, (g^, q^) — > {q^,Pb, u^), to write 
Pa{q,q,t), Pm^ fm(q,Pb{q,q,t),t^ , u'^{q,q,t). 

Provided with these functions, pA{q, q, t), we define 



l{q, q, t) {q^PA - H(p, q, t) 
to recover the Lagrangian dynamics. 



V 



q^PA - HTip,q,u,tj^ 



V 



(3.38) 



(3.39) 



dL{q, q,t) 
dq^ 

dL{q, q,t) 
dq^ 



^ f r_^MBqBdpB _ dpBdHr _ 

\ dq"^ dpB dq^ ) 



PA{q,q,t) , 



(3.40) 



dHT{p,q,u,t) 



V 



dq' 



V 



(3.41) 



along with the M constraints, Pm = fm{q^:Pb: t) ■ 



This analysis shows the equivalence between the Hamiltonian and the Lagrangian formalism, 
up to the technical assumption on the existence of the inverse map (3.38).^^ 



Notice that the equivalence between the Lagrangian and total Hamiltonian system can be shown in a large 
number of ways, see e.g. the exercises 1 .2 and 1 .4 of the book [7] for some alternatives. The latter exercise is based 
on the general result about the elimination of "auxiliary fields" (in the present case, the 'auxiliary' fields are pb and 
u"*) via their own equations of motion . 
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4 Total Hamiltonian dynamics 
4.1 Poisson bracket 

On the 2A/^-dimensional phase space, with coordinates (q'^tPb), equipped with the Z2-grading, 
we define the Poisson Bracket as 

^ '^^^-^^ ^ ^ dq^dpA ^ ^ dpAdq^- ^ - ^ 

The Poisson bracket can be rewritten in a more compact form, 

\ I 51 - (4 2) 

where the arrows indicate the direction the derivatives act. It satisfies the graded skew-symmetry 
property, 

[F, CIp.^. = -(-!)#-#- [G, F}^^ , (4.3) 

the Leibniz rule, 

^^^}p.B. = [F, G^}p.i.. K + {-l)*-*-G [F, K}^^^^ , 

(4.4) 

[FG, K}^^ = F [G, K}^^ + {-l)*o#- [F, K}^^ G , 
and the Jacobi identity, 

[[F, G}p^. , K}^^ = [F, [G, i^}p.5.}p.^. - [G, [F, K}^^ }^^ , (4.5) 

or equivalently, 

[F, [G, i^}}p^. + [G, [H, F}}p^. + [i/, [F, G}}p^. = . 

(4.6) 



Let 1^ denote the Hermitian conjugation such that (a by = 6^ , /.e. ^ is an involution on the 
algebra of functions on the phase space. Reality condition on the phase space reads, 

g^t^g^, p^^^{-l)*Ap^, (4.7) 

because the symplectic form q'^pA must be real and 

{q^'pAY = p\{qy = i-lf^iq^pl ■ (4.8) 
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Hence we have 



and hence. 



(^] ' - ^ ' - ^ (4 9) 

\dq^J ^ ' dq^ ' {dpAj ^ ' dpA ' 

{[F,G}p^^y = [Ft, Gt}p.s. = -[G\F^}p.B.. (4.10) 



4.2 Time derivatives - preliminary 

With a given total Hamiltonian, 

M 



HT{p,q,u{p,q,t),t) ^ H{p,q,t) + J2 <l>m{p,q,t)u'^{p,q,t) , (4.11) 

m=l 

and the generalization motivated at the end of Subsection 3.4, we let the following formulae 
govern the dynamics of the full phase space, 

q^ = {-l)*-—L, p^ = --4, (4.12) 

OPA oq^ 

where the variables u^{p, g, t) are to be understood as functions of p, g, t, of which the explicit 

forms are not yet specified. Note that with the restriction on the hypersurface V, i.e. putting 

(pra = after taking the derivatives, the above equations reduce to (3.37,3-40,3.41) which already 

indicates some equivalence between the total Hamiltonian dynamics and the Lagrangian dynam- 
ic 

ICS. 

The time derivative of an arbitrary quantity F{p,q, t) takes a simple form^ in terms of the 
Poisson bracket (4.1), 

dt -[^'^^>P.B. + ^- (4-14) 



^^It is worthwhile to note that the above equations (4.12) indeed govern the full dynamics of all the coefficients 
[pA\j and [q^]j of the Grassmann algebra (see Appendix B). 

^"^One may generaUze Eq.(4.14) by adding an arbitrary quantity proportional to the constraints to the right hand 
side, 

— -[F,Ht}p.^. + — 

See (4.48) and also the Dirac bracket (6.10) for further discussion. The Dirac bracket gives an alternative dynamics 
which coincides with the dynamics of the Poisson bracket only on V, but it is more suitable for quantization. 



[F, ^tIp.b. + — + ^mC^iF) ■ (4.13) 
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The crucial viewpoint we adopt here is the following: Though the above equations supple- 
mented by the primary constraints are equivalent to the Lagrangian formalism (as shown in 
Subsection 3 .4) we consider them to be more fundamental. Namely, without restriction on the 
hypersurface, we let them govern the entire phase space. Then, we try to make sure that the dy- 
namics can be consistently truncated to the hypersurface. Namely we will look for the necessary 
and sufficient conditions to maintain the constraints throughout the time evolution. By imposing 
so, it may well be the case that we determine some of the unknown variables m"* completely, 
at least the values on the hypersurface, and obtain further consistency conditions or 'secondary 
constraints'. In the latter case, both the primary constraints and the secondary constraints should 
be imposed to define the hypersurface, say V. In this way, we indeed make our Hamiltonian dy- 
namics be consistent with the Lagrangian dynamics if one considers the space of physical states 
to live on the constraint surface V . 



4.3 Preserving the constraints - primary and secondary constraints 

On the hypersurface V, we have 



[F,H}p_^+[F,cf>^} 



P.B. 



OF 
'dt 



(4.15) 



where we introduced the notation ~ for the "weak equality" defined via the equivalence relation 
f^g ^ f\v^g\v ^ f^9 + <f>mz'^. (4.16) 

We remind the reader that \v means the restriction on the primary constraint surface V, or the 
expression of (p, q) in terms of the variables as in (3.16), after taking the partial derivatives. 

In other words, the symbol ~ stands for "equal on the primary constraint surface". This symbol 
makes the distinction with the "strong equality", which is the usual equality throughout all phase 
space, and proves to be convenient in order to avoid repetitions of the restriction on V for every 
term in a lengthy expression. On the left hand side, we get from (3.18) that 



V 



^( F{p,q,t)\y + (l>mF''' 



V 



^(F(p,?,i)y + 0^F^ 



(4.17) 



For the consistency with the Lagrangian formalism, the time derivative of the primary con- 
straints must vanish on V, 



dt 



(^[<Pm, H}p^ + ^ + I*^-' '^n}p.s. - ' (m, n = 1, 2, • • ■ , A<) , (4.18) 
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which precisely corresponds to (3.29). We focus on the following supermatrix defined on (q"^, t). 



^^mn:= [0m,0n}p.B.|y = -(-1)"" [</'n>m}p.B.|y= ( ) , (4-19) 

where A± are symmetric or anti-symmetric (A^ = ±A±) bosonic matrices while ^ is a fermionic 
matrix. Now Eq.(4.18) can be taken as a set of linear equations with the M unknown variables 
u'"(l<m< M), 

[0„, H}p^ + ^ - -^nmU"^ , (4.20) 

where, surely, the lefthand-side and il are given by fixed functions on V of either (x, t) or, equiv- 
alently, {pb, q^,t)P 

Let us analyze the process more concretely. Without loss of generality, assuming that all the 
primary constraints are real, 

= (l>^ , [u^V = {-l)*"^u^ . (4.21) 

the supermatrix is anti-Hermitian, Vtmi = — (i7;m)^ so that all the matrices are real, A± = 
A'^, = The forthcoming analysis is rather technical and complicate since one includes 
fermions. The main results are summarized in Subsection 4.4 to which the reader may jump di- 
rectly in a first reading. 

Under the real linear transformation,^^ 

{L*)k"' = (-l)#-(#fc+#-)(Lfc"*)* = Lk"" , 

<t>m Ljct>k = ct>k{L^fm, (4.22) 

, ((L^)-l)"^fc^i^= (-l)#-+#*-f/'=(L-l)fc"^, 

the contraction is invariant, the reality condition (4.21) is preserved, and the supermatrix 

transforms as 

[4>A}p.b\v ^[</'»p.B.ly^^> (4.23) 



^^For a given set of the generating elements of the Grassmann algebra A ^ (as in Appendix B) the above formula 
(4.20) can be, in principle, completely analyzed. Furthermore, for the consistency of the Lagrangian mechanics, there 
must be a solution thereof. Indeed, we always impUcitly assumed that we disregard any inconsistent Lagrangian Uke 
L = q which would lead to = 1 = 0. In the same spirit, we can also expect that, if necessary, there may occur 

some more extra constraints on V , that is 'secondary constraints'. 

26Note that, in general, = (-l)#".(#-+#fc)Lfc'», (L*)fc'" = seeEq.(C.2). 



27 



where we have set L\y — L for simplicity. As shown in (D.31), one can always transform any 
anti-Hermitian supermatrix into the following 'canonical' form by a real linear transformation,^^ 



(4.24) 





( 











\ 


LVtL^ = 





s_ 





V' 

















^ 











where all the matrices are real, b± — 




= 4> 




-- r ; 


h 



matrices so that h± ^ exist ; s± are bosonic products of fermions {i.e. even number of products 
of fermions); and h± — ±6^, s± — . It may be the case that s± and/or vanish. 

With the decomposition of the index m into c, c for bosonic variables and a, d for fermionic 
ones, which should be obvious from the inspection of Eqs.(4.25)-(4.28), the consistency condition 
(4.18) now splits into 



+ 



and 



dt 
dt 

d(j)c 

dt 



(4.25) 
(4.26) 

(4.27) 
(4.28) 



The meaning of the first two equations, (4.25) and (4.26), is clear. Since b± are nondegener- 
ate, they fix the unknown variables, u^, u", completely as functions of {q^,Pb, t) or {x, t) on the 
hypersurface V . 



The analysis of the last two equations, (4.27) and (4.28), is somewhat tricky. Before the full 
analysis, we first focus on the purely bosonic systems, which was the case studied by Dirac [4]. 

^^Contrary to the usual complex number valued Hermitian matrix, a Hermitian supermatrix may not be completely 
diagonalizable. However, if the Hermitian supermatrix is nondegenerate, it is diagonalizable. See our Lemma 4 in 
(D.31). 
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Bosonic systems. 

In the bosonic systems, the equations (4.26) and (4.28) simply do not appear, and the es- 
sential relations are 



ii ^}p.B. + 



at 










[ 1 



(4.29) 



Since 5_ is nondegenerate, the variables if are completely determined on V in terms of 
the variables (pa, Q^^t), while the other variables remain as locally free variables at this 
stage. The vanishing of the second row in the left-hand-side can give some new, namely 
'secondary', constraints. In this case, the primary and secondary constraints define together 
a smaller hypersurface, say V C V, and some of the p^'s can be expressed in terms of 
others. Then the already determined variables u'^ should be further restricted on V', making 
the first row hold on V too. We note that the number of secondary constraints, say n, are 
not greater that the number of the yet free variables, m = dim{M' }, n <m. 

The next step is to consider the time derivatives of the secondary constraints, analogously 
to (4.18). Regarding them as linear equations in the variables u'^ and taking some linear 
transformations to the canonical form, (D.21), 





H 











)l 





< A; < n < m , 



0(n— fc)xfe ^{n—k)x{m—k) 

(4.30) 

one can determine the variables u'^' completely on V, and there may appear new, namely 
"tertiary", constraints x'\y, ^ 0. Again the number of the tertiary constraints, are not 
greater that the number of the surviving locally free variables u'^' , since (n — k) < (m — k). 
The procedure may go on, but it should terminate at certain point, since the total number 
of constraints should not exceed the dimension of the whole phase space for any consistent 
dynamics. By a slight abuse of terminology, one refers to all these new constraints as 
'secondary'. 

Eventually, we end up with a set of constraints, 

[(ph^O,l<h<M + M'}, (4.31) 

of which the ranges 1 < h < M. and M. + 1 < h < M. + M.' respectively correspond to 
the primary and secondary constraints. They define the hypersurface, V, 

V := [ip,q) \ (ph = , 1 < h < M + M' Y (4.32) 
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All the constraints 0/i (1 < /i < + M.') are static on V, in the sense that the restriction 



on V is preserved by the time evolution, 



. In other words, they satisfy 



- H}p.B. + ^) = E ( b/., Mp.b. «™) 

^ /V m=l 



X 




(4.33) 



in the canonical form on V . 



Some of the variables, i.e. the m™' 's, are completely determined^^ on V as functions of 
{p, q, t), while the others, the u"" 's, if any, remain as locally free (that is, arbitrarily time 
dependent variables). The latter correspond to the zero eigenvectors of the {M.+M.') x M. 
matrix, [4>h,4>m}p,B.\v ■ 



Generic systems. 

Now we return to the equations, (4.27) and (4.28), which are relevant to the generic systems 
of bosons and fermions. 



U TJX J. ^"f"^ 
y [0a,-H|p.B. + J 



is-) 



ab 



^a0 



(4.34) 



In general, the complete analysis of the above formulae is always possible, if we introduce 
explicitly the Grassmann algebra basis of (B.2). By expanding all the quantities in terms of 
the basis accompanied with the real or complex number coefficients, i.e. [q^]j, \pa]j, one 
can convert them into the linear equations in [u°-] j, [u°'] j, over M or C. The linear equations 
can be completely analyzed, essentially in the same way as in the previous bosonic case. 
After all the finitely repeated procedures, the results will be parallel : some of the coeffi- 
cients, are completely determined in terms of \pa]j, t), while others 
remain as locally free parameters, implying that the time evolution in the phase space is not 



^**Strictly speaking, what we have determined are the variables on the hypersurfaces, {m™ }. For the generic 

dynamics in the full phase space, we may either employ them literally as they are, or use the continuously extended 
functions which have nontrivial dependence on the orthogonal directions to the hypersurfaces. In any case the 
dynamics on V is the same. 
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deterministic. There may appear secondary constraints in terms of the coefficients, [?^] j, 
\Pa]j. 



However, in practice we favor the Lagrangian systems which do not require any explicit 
use of the basis for the Grassmann algebra. In such 'good' systems of both bosons and 
fermions, all the expressions can be written collectively in terms of (p, q, u, t), rather than 
([p] J, [q'] J, [li] J, i), as if in the bosonic system. We summarize the results in the following 
separate subsection. 



4.4 Hamiltonian dynamics after analyzing the constraints - summary 

In all the bosonic systems as well as all the 'good' systems for bosons and fermions, in the sense 
that the explicit use of the basis of the Grassmann algebra is not required, we have the following 
generic situation: 

• The Hamiltonian, H{p,q,t), as well as the primary constraints, (f)„i{p, q,t), are given as 
functions on the 2A/'-dimensional full phase space with coordinates {pa, q^) with 1 < 
A,B<M . 

In particular, the primary constraints define a {2J\f — Al) -dimensional hypersurface 

V ^{{p,q)\(i)m{p.q,t) = Q, l<m<M). (4.35) 

• The total Hamiltonian is a sum of the canonical Hamiltonian and linear combinations of 
the primary constraints, 

M 

Ht{p, q, u{p, q, t),t) = H{p, q,t) + J2 Mp, Q, t) u'^ip, q, t) . (4.36) 

m=l 

• The dynamics of the whole 2A/'-dimensional phase space is subject to 

^' = i-lf^f^, PA--^, (4.37) 

dpA oq^ 

so that the time derivative of an arbitrary function, say F{p,q,t), on the 2jV-dimensional 
phase space reads 

dF(p,g,t) _ , dF 

where the variables u^{p, q, t) are not yet specified, but by looking for the necessary and 
sufficient conditions to maintain the hypersurface V , in order to be consistent with the dy- 
namics, we may determine some of them completely. In general we encounter the following 
situation: 
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• The whole primary constraint surface V may not be consistent with the dynamics, in the 
sense that it may not be preserved by the time evolution. It may well be the case that only a 
subset of V, say V C V,is preserved. The constraint surface V is specified by the primary 
as well as the secondary constraints, 

V = {(p, q) I Mp, q,t)~0, l<h<M + M'}, (4.39) 

where {<ph} denote the complete set of constraints indexed by h , such that A4 + 1 < h < 
M. + M.' correspond to the secondary constraints. We introduced the notation Ki for the 
other "weak equality", that is defined as 'equal on the constraint surface V, 

f^g ^ f\v = 9\v ^ f = g + <Phz''. (4.40) 

We emphasize the important distinction between the two weak equalities, because (4.16) 
implies (4.40) but the converse is not always true since V C V . All the constraints, in 
principle, can be solved to express any point on V by 2J\f —M. — M' independent variables, 
say y' (with l<i<2Af-M- M'), 

V = {{p,q) = {f{y,t),g{y,t))} , (4.41) 

which provide a local coordinate chart on V. On the other hand, A4 + M.' of all (f)hS can 
be taken as complementary coordinates for the orthogonal directions to V in the full phase 
space. The entire 2jV-dimensional phase space has then two sets of coordinate charts, 

{pA,q^) , l<A,B<Af 

(4.42) 

<^ {y\ (f)h) , l<i<2J\f-M-M', l<h<M + M'. 



• There exists at least one set of solutions for { m™, 1 < m < M.} and a {M. + M.') x {M. + 
M.') supermatrix satisfying for all the constraints, 

= [<l>h, H + (Pmu'^yp,^. + ^ = <l>h'T''\ , l<h<M + M', (4.43) 
or equivalently 

H}p^ + ^ ~ - t'^'^' ■ (4-44) 

where indicates that the equality strictly holds after the restriction on V or, equivalently, 
putting {p, q) — {f{y, t) , g{y, t)), after taking the derivatives. 
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The most general solution of Eq.(4.44) reads 

u^ip, q, t) = U"^{p, q, t) + V"\ip, q, t) v\t) , (4.45) 

where (with I < i < T < A4) are arbitrary time dependent functions, U"^{p, q, t) is 
one particular solution, and V^i{p, q, t) span a basis of the kernel of the [M + M') x M 
supermatrix [cph, (pm] p,b. I v ' 

M 

XI [ <^m}p.B. ^""i ~ , for all 1 < /i + . (4.46) 

m=l 

• When (pm—Pm — fm{Q: Pa: t)-, Substituting the most general solution (4.45) for into the 
expression of the velocity 

q^ = q^{q'',Pa,{-lf-u^{p,q,t),t) 
given by (3.27), the total Hamiltonian reads simply: 

Ht{pa, q"", t) = q-'pA- Hq"" , q^, t) . (4.47) 

• As in (4.13), one can generalize the total Hamiltonian dynamics by adding terms propor- 
tional to the constraints. In particular, still preserving the Poisson bracket structure of the 
time evolution, (4.38), - which is essential for the quantization - one can modify the to- 
tal Hamiltonian alone by adding freely terms quadratic (or higher) in the constraints, both 
primary and secondary, 

HT^H + (l>mU"' =^ HT^H + (l>mU"' + lMh'w''''' . (4.48) 

where 1 < h,h' < M + M', while 1 < m < M, and w'^^' = (^—l)#h*h'w'^''^ are newly 
introduced local parameters, being arbitrary time dependent functions. The modification 
does not affect our previous analysis at all, and the resulting dynamics remains the same on 
the hypersurface V . 

• Other characteristic features of the total Hamiltonian dynamics are discussed in the follow- 
ing subsections. 
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4.5 First-class and second-class 

We define a dynamical quantity, say g, t) a function on phase space, to be "first-class", if it 
has zero Poisson brackets with all the constraints, both primary and secondary, on V, 

F,(j)h} ^0, l<h<M + M'. (4.49) 

L J P.B. 

Otherwise it is said to be "second-class". Physically, this distinction is extremely important, par- 
ticularly for the constraints, as will be explained later on in this subsection. 



4.5.1 Main properties of first-class quantities 

We enunciate some useful properties of first-class quantities: 

• If the function F on phase space is first-class, then the Poisson bracket [F, (ph}pB ^ust be 
a linear combination of the constraints, 



""Ah] 



P.B. 



fh^'cl>h'. 



(4.50) 



IfF is first-class, then for any arbitrary dynamical variable, say G{p, q, t). 



'f.G\v] 



P.B. IV 



P.B. IV 



(4.51) 



Proof: For an arbitrary analytic function, G{p,q,t), with the new coordinates for the 2J\f- 
dimensional phase space, (y*, (ph) from (4.42), we define G{y,(f),t) := G{p,q,t). 
Then from the property of being first-class and the Leibniz rule (4.4) of the Poisson 
bracket, one can derive the result, 

L J P.B. 

(4.52) 

□ 



F,Giy,<P,t) 



P.B. 



F,G{y,0,t) 



F,G\ 



P.B. 



P.B. 



• The Poisson bracket of two first-class quantities is also first-class. 
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Proof: This can be shown from the Jacobi identity (4.6). For two first-class quantities, 
say Fi and F2, 

J P.B. V 

J P.B. V L J P.} 



Fi 



P.B. 



JP.B. Iv 



(4.53) 



□ 



4.5.2 First-class constraints as gauge symmetry generators 

From Eq.(4.46), the complete set of primary first-class constraints is given by 

In virtue of (4.48) and (4.54) the total Hamiltonian reads now 
where H' is defined as 

H':^H + 4>mU'^ , 
and satisfies, for any 1 < h < A4 + A4', 



d(ph 



P.B. 



dt 



(4.54) 

(4.55) 
(4.56) 

(4.57) 



Thus, when the constraints do not have any explicit time dependence, both H' and Ht are first- 
class, and up to quadratic terms in the constraints the total Hamiltonian is a sum of the first-class 
Hamiltonian H' plus a linear combination of the primary, first-class constraints with arbitrary 
functions of time as coefficients. Notice that such a decomposition is not unique since U"^ can be 
any solution of the inhomogeneous equation (4.44). 
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The time derivative of an arbitrary function q, t), c.f. (4.38), on the 2j\/'-dimensional 
phase space can be rewritten as 



dF{p,q,t) 
dt 



P.B. 



V\t) + 



OF 
'dt' 



(4.58) 



As a result of the appearance of the arbitrary time dependent functions v'^{t) , the dynamical vari- 
ables at future times are not completely determined by the initial dynamical variables. 

We should recall the following viewpoint by Dirac [4]: "all those values for the p's and g's at a 
certain time which can evolve from one initial state must correspond to the same physical state at 
that time". A natural definition of the space of physical states is thus as the set of initial variables 
that should be given at some given moment of time, say to, in order to determine completely 
the time evolution via the equations of motion. The presence of arbitrary functions of time, 
v'^{t) , in the total Hamiltonian Ht signals that the phase space contains some unphysical degrees 
of freedom. Indeed two different choices of arbitrary functions, say v\ and V2 , would lead to 
distinct total Hamiltonians and thus different time change of a dynamical variable, say F. After 
some time interval St , the evolutions of F would differ by 



The key philosophy we stick to is the standard one that different choices of the local gauge pa- 
rameters correspond to different total Hamiltonian systems, which nevertheless should be taken 
equivalent, describing the same physics. Following the viewpoint advocated in subsection 2.2.2 
for Lagrangian systems, this means that Eq.(4.59) defines an ambiguity in the time evolution that 
should be physically irrelevant. In other words, the transformation (4.59) is a gauge symmetry. 
In modem terminology, this implies that: 

(i) a physical state is represented by an equivalence class, where one mods out by the gauge 
symmetries, therefore 

(ii) the space of physical states must be understood as the quotient of the constraint surface by 
the gauge orbits and 

(iii) an observable is a gauge invariant function on the constraint surface. 

The space of physical states is also a symplectic manifold^^ and is sometimes called "reduced 
phase space." As a short dictionary for physicists on mathematical jargon, we may say that the 

■^^The section 1.4.2 of [7] is devoted to the subtle counting of physical degrees of freedom in the Hamiltonian 
context, which is equal to half the dimension of the symplectic manifold. 



SF= [F,cl>]^'}ivl-v\)St. 



(4.59) 



36 



former quotient space is usually called "symplectic reduction" by mathematicians while they 
would introduce gauge transformations via a "Lie group action" and might refer to Noether's 
theorem as the "moment map" (see e.g. [13] ). 

A remarkable property is that the Poisson bracket {H', (f)}^^} between the first class Hamilto- 
nian H' and any primary first class constraint (j)}^^ is also a gauge symmetry generator. This can 
be shown by comparing the time evolution successively determined by (i) the total Hamiltonian 
Ht during an interval Sti and after by the first-class Hamiltonian H' during an interval 5t2, or (ii) 
the same operations, but in the reverse order. The net difference must be a gauge transformation 
since Ht and H' define the same evolution of physical states. By using the Jacobi identity, one 
may check explicitly that this gauge transformation is indeed generated hy {H', (p]^^} . 

As one can see in (4.59), the primary first-class constraints generate gauge symmetries. A 
natural question is whether the converse is true: are all gauge symmetries generated by primary 
first-class constraints? In full generality, the answer is no. This can be understood from the fact 
that the Poisson bracket {H', (p}^^} is also a gauge symmetry generator. From Eq.(4.50) we know 
that this Poisson bracket is a linear combination of first-class constraints, but it is not guaranteed 
that only primary constraints appear. Therefore, some secondary constraints may also generate 
gauge transformations. Then another question arises: do all the secondary first-class constraints 
generate gauge symmetries? Dirac conjectured that the answer would be yes. But, again, in full 
generality the answer is no, although in most physical applications the answer is yes.^'' This is 
the reason why first-class quantities have such a distinct status. 

4.6 Extended Hamiltonian dynamics 

One can generalize the total Hamiltonian system further to the so-called "extended Hamiltonian" 
system. We define the extended Hamiltonian in a similar way to the total one except that the 
former includes all first-class constraints (the primary as well as the secondary, tertiary, etc), 

He{p, g, t, v\ v^\w) := Ht + 0f = H'{p, g, t) + cj>l^'v%t) + cj>l;^V (t) + ^Mh'W^"' , 

(4.60) 

where (with 1 <i' <T < M') denotes the secondary first-class constraints. 

The extended Hamiltonian is usually preferred because, from the Hamiltonian point of view, 
the distinction between primary and secondary constraints is actually irrelevant. The distinction 
becomes important only if one wants to make contact with the Lagrangian formulation, as in 

'"'a counterexample of the Dirac conjecture is given in subsection 1.2.2 of the book [7]. A proof of the Dirac 
conjecture under some hypotheses is provided in its subsection 3.3.2. 
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Subsection 3.5. So, one may let He governs the whole dynamics rather than Ht, 

dF dF 

-^ = [F,H,}...^-. (4.61) 

Compared to the total Hamiltonian dynamics, the constraint surface V, is still preserved but a 
generic other object, say F (rather than the constraint (ph), undergoes a different time evolution 
from the total Hamiltonian dynamics, even on V. Indeed, [F, 4>\f}p.B.\^^; 7^ in general. Thus, 
unlike the total Hamiltonian dynamics, the extended Hamiltonian dynamics is in general different 
from the Lagrangian dynamics. Still, if the Poisson bracket of a quantity F with any secondary 
first-class constraint is zero, then its evolution on V are the same. 

With an arbitrary local parameter e{t) , if there exists a gauge symmetry generator Q}^^ — 
Q^^^ e{t) which takes any solution {q^,PB){t) of the extended Hamiltonian dynamics to another 
by 

V = [ Q'''}p.B. e{t) , 5pB^[pB, Q'''}p.B. e{t) , (4.62) 
then at any time, say to, one can start to transform the solution to another without changing the ini- 
tial data {q^:PB){to) by simply setting s{to) = 0. Thus, again the future dynamical variables are 
not uniquely determined by the initial data. From the point of view of the gauge symmetries, the 
main difference between the total and extended Hamiltonian dynamics is that the latter assumes 
the Dirac conjecture applies. In this case, a dynamical quantity F{q^,pB, t) on the phase space 
defines an observable if and only if its Poisson bracket with any first-class constraint vanishes 
weakly 

[ F, (f)]f^] ^Q, {l<i" <I + I'). (4.63) 

Notice that an observable has been defined as a function on the constraint surface, so one should 
identify two functions that coincide on V , i.e. the observable corresponding to the first-class 
quantity F is the equivalence class for the weak equality. The conclusion is that the physical 
quantities (that is, the observables) undergo the same evolution under the total and extended 
Hamiltonian dynamics.^^ 



4.7 Time independence of the Poisson bracket 

For a given set {v^{t)} of the local functions, the dynamical variables (p, q) follow a unique and 
invertible trajectory in the phase space such that there exists a one to one map between (p, q) and 
the initial data, 

PA(^,Po,go), pa(^o,po,?o) =poa; q^(t,po,qo), q^{to,Po,qo) ^qi ■ (4.64) 

^'For more comments on the relationship between the total and extended Hamiltonian formahsms the reader may 
look, e.g. at [14]. 
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A crucial fact follows, proven in Eq.(A.15): The Poisson bracket is independent of time. 



L ' iP.B. ~ L ' iP.B.ho 

or 



(4.65) 



-1 ^ ^ = -1 1, 7^— 7^ • (4.66) 

dq^ dp A op A oq^ dq^ dpoA opoA dq^ 

Namely the time evolution generated by the total Hamiltonian is a symplectic transformation. 



4.8 Other remarks on the total Hamiltonian formalism 



A useful identity. 

For an arbitrary function F{p, q, t) , we have the following identity,^^ 

d 



(4.67) 



proven in Eq.(A.13). In words, the following two actions, taking the time derivative and 
restricting on V, commute with each other. Intuitively, this is obvious since we imposed 
that V be preserved under the time evolution. 

The time derivatives of the primary first-class constraints (pj^^ are of the general form, using 
Eqs. (4.43) and (4.55), 



1st 



dt 



•If - } 



P.B. 



v'v^^^ = , 



(4.68) 



P.B. 



for all h^l,2,---,M + M' . 



Namely the time derivative of any primary, first-class constraint decomposes into two parts, 
one being independent of the local gauge parameters v'^{t) , and the other one being first- 
class and orthogonal to the local gauge parameter. 



^^The above relation (4.67) should be compared with 

^ dpA f dF 
y dt \dpA 



d_ 
dt 



dF 
'dt 



+ 



dt \ dq 



dF 
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• In the static case, where there is no explicit time dependence in q) and H'{p, q) de- 
fined in (4.56), the time derivative of any primary as well as secondary first-class constraint, 
0^***, is first-class^^ too, since 0^''* = [0^''*, Hrjp^ so that 

^4 = [01^*, [Ht, 04}p.B. - [Ht, [<i>''\ 04}p.B. ^ • (4-69) 

J P.B. 

Hence , in the static case, for all the first-class constraints, both the primary ones, and 
secondary ones, (t)\f, we can write 

(4.70) 

• Combining the above two results in the static case, [0^^®*, H'}p_b. and v^f are separately 
first-class constraints. 



^Actually, this is an example of the fact that the Poisson bracket of two first-class quantities is also first-class 
(4.53). 
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5 Symmetry in the total Hamiltonian system 
5.1 Symmetry from the Lagrangian system - revisited 

In the Lagrangian formalism,^'^ the notion of symmetry corresponds to a change of variables^^ 
q'^ which leaves the Lagrangian invariant up to the total derivative term as in (2.14), 

<\K 

L{q',q',t)=L{q,q,t) + — . (5.1) 

As a consequence, the symmetry takes one solution of the Euler-Lagrange equations to a new one 

as in (2.21), 



dL{q,q,t) ^ d_ r dL{q,q,t) \ ^ dL{q',q',t) ^ d_ r dL{q',q',t) 
dq^ dt \ dq^ J dq'^ dt \ dq'^ 



(5.2) 



As discussed in Section 2.2, if the infinitesimal symmetry transformation ^ + 5q'^{q, q, t) 
depends on g"^, g"^, t only, then the quantity 6K in 5L = ^SK must also depend only on q^, q^, 
t , and hence so is the corresponding Noether charge: 

Q{q, q, t) = Sq^{q, q, t)pA{q, q, t) - 6K{q, q, t) . (5.3) 
5.1.1 Change of variables 

Henceforth in the present subsection, i.e. until Eq.(5.15), we take (g, pa, g"^, t) as the independent 
variables for any quantity which carries a hat symbol. For instance, we set 



Pa — ( ^%^) . - „ . > Pa-^Pa, Pm fm{q^ , Pa, t) , 

Tq^{q^,Pa,cr,t) ■.^5q^(q^,h%q,pa,r,t),q'^,t) , (5.4) 

mq^,Pa,q'^,t) := SK^q^ ,h%q,Pa,q\t),q^,t^ , 

where we substituted g" by h'^{q^,Pa: t) according to the relation in Eq.(3.7). Notice that 

dpB 



dq''- 



0, (5.5) 



^''This subsection is parallel to Sec. 3.4 where the Hamiltonian corresponds to the time translational synametry 
generator. 

^^In order to avoid confusion with some subsequent notations, we slightly changed the convention by using 'prime' 
instead of 'tUde' to denote the new variables. 
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as follows from (3.7). Furthermore, in agreement with Eq.(2.38) we define another function 
depending on (f", t). 



SpA{q,Pa,t,q" 



dL 



d{Sq^) fdfdL 



dL 



dq^ \dt \dq^ J dq^ 



(5.6) 



d{6K) d{Sq^) dL 



dq^ 



QqA QqB 



which has been defined in such a way that the r.h.s. is independent of the accelerations, as it 
should in the Hamiltonian formalism to come. Similarly, notice that the evaluation of (2.34) at 
= i), gives 



( d{SK) \ ^ f d{6q^) 
V dq^ ),^=H^ V dq' 



A I Pb, 

q°-=h°- 



(5.7) 



From the next equation (5.10) until the equation (5.15), the partial derivatives acting on any 
quantity with a hat symbol, say 



A{q^,Pa,cr,t):^(A{qA,t)) . , 

\ / q''=h°-(q,pa,q'^,t) 



(5.8) 



is taken regarding {q^^ps, (f^, t) as independent variables, while for unhatted quantities the in- 
dependent variables are (g^, q^ ,t). Concretely, this means that one should make use of the chain 
rule, so that 



dA 
dq^ 

dA 
dq^ 



dA\ dh^ /dA\ 

dq'^Jq^=h^ dq^ \d^^^ q^=h^^ 
dh^ /dA\ ( dA 



dq^ V5g&/ga=/ia 



+ ^ 



\ dq^ J Q^=h^ 



(5.9) 



On-shell, the variation (5.6) is equal to the variation of the momenta. Moreover, it satisfies 

d SpA{q, Pa, cT, t) d 5q^{q, Pa, cf^, t) ( d^L{q, q, t) ' 



dq'^ 



d(f' 



QqBQqA 



q^ = h^{q,Pa,q'^,t) 



(5.10) 



dq 



dq' 
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The proof of these two equalities is provided in Appendix A. Similarly to the corresponding 
steps, we also have the following identities: 

f_l)#A#S "1 ^ Q 

^ ^ dq^ dpAdh^dq^ 

(5.11) 

PC dq"^ dh'^dq^ dq^ dh"- dpc 



The above substitution induces a Noether charge depending on {q,pa,t): 

Qiq^,Pa,t)^Tq^PA-'5K. (5.12) 

From (5.7) one can easily verify that the Noether charge is indeed a function of q"^, Pa and t only 
(i.e. it is independent of q^), 

dQ ( dh' d{6q^) ^ d{Sq^) \ ^ fdh^ d{5K) ^ d{5K) \ ^ ^ . ^ 



dq"^ ydq"" dq"- dq"" J ydq"" dq"- dq"" J 

Furthermore, with (5.7) we get 

(5.14) 

= (-!)#« + (-l)#a#m'Sqmdfrn ^ 

OPa 

and with (5.6), 

dQ{q,Pa,t) _ r d{6q^) ^ dh^ d{dq^) \^^ ^ djdK) dh^ d{6K) 

QqA y QqA q^a gj^a J q^a q^a q^a qj^o. 

^-TpA + {-l)*-*-Tq"^^. 

In order for the Noether charge Q{q,Pa,t) to generate the symmetry transformations {5q, 5p) 
via the Poisson bracket in the corresponding Hamiltonian system, the derivatives of in the 
r.h.s. of (5.14) and (5.15) should be absent. Instead, in a spirit similar to the total Hamiltonian, 
we will define a "total Noether charge" which will be a function of (p^, t) rather than^^ 



36Note that the equation (3.27) impUes that ^ = = (-l)**^^^*'"^!^^!? 
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5.1.2 Total Noether charge 

Let us denote by v"^ the explicit expression of q"^ in terms of the variables , Pa, q"^ and t , as 
in (3.27). Substituting the velocities q by their explicit form v{q,Pa, q^, t) henceforth, we take 
{q^,Pa, Q^, t) as the independent variables for any quantity which carries a tilde symbol. In other 
words, we set 

PA-^(^%T^)^ ^^PA, Pa-=Pa, Pm-=fm{q,Pa,t), 

jKiq, Pa, t) := SK{q, v, t) = lK{q, Pa, v^, t) , 

Sq^iq,Pa,q"',t) -.^ 5q^{q,v,t) ^'5q^{q,pa,v'^,t) , 

SpA{q,Pa,q"',t) := SpA{q,Pa,v'^,t) . 
The total Noether charge is then defined as 

QT{p,q,q'^,t) := Sq^{q,Pa,q'^,t)pA- SK{q,Pa,q"',t) 

= Qiq,Pa,t)+ 3q''{q,Pa,q"',t)(^Pn- fniq,Pa,t)Y 



(5.16) 



(5.17) 



Contrary to the quantity Q{q,Pa: t) in (5.12), in the above expression of Qt(p, q, t, (f) the con- 
strained momenta Pm have not been substituted by the primary constraints pm — fm{q,Pa, t) , 
since it is the untilded momenta which multiplies Sq . It follows from Eqs.(5.13)-(5.15) that^^ 



dQTip,q,t,q"') ^ d{6q'') 



(Pn- fn{q,Pa,tfj , 



dQT(p,q,t,q"') .-si^ATA ^^iSq"") f \ , j.. 

= {-l)^^6q +— [Pn- fn{q,Pa,t)] , (5-lo) 

dpA opA V / 

dQTip,q,t,q'^) ~ d{6q 



dq"^ dq 



~ d[ dq ) ( \ 

- SpA + —g-A~ [P'' ~ fn{q,Pa, t) j . 
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'Similar equations to (5.18) can be straightforwardly obtained either for the case where we take {q,Pa, t, g™) as 
independent variables or for the case where all the = (— are completely determined in terms of {q,Pa, t) 
and free parameters v^{t) after solving all the constraints (4.45). 
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In particular, in terms of the Poisson bracket we have 

[QT,q^}p.B=-Tq'' 

[Qt-,Pa}p.b. = -5pA 

Note that the expressions are consistent with the integrability relations e.g. dp^dqmQx — 
{—l)*^*"^dqmdp^QT, thanks to (5.10) and (5.11). Obviously on the primary constraint surface 
V, where pm — fm{Q,Pa, t) , the above relations get simplified: the first relation in (5.18) means 
that Qt becomes independent of q"^ onV , and the other equations lead to 

[ Qt, /}p.s. ^ -Tq^ , [Qt, Pa}rb. ^ -TpA . (5.20) 

In the Hamiltonian formalism to come, the relations (5.20) will be interpreted as the property 
that the total Noether charge Qt generates the infinitesimal symmetry transformations on V and 
on-shell.^^ 



[g^, dq'^jpB.^Pm - fm{q,Pa,t)^ , 
[PA, Sq"^}p.B. (pm - fm{q,Pa, ^)) ■ 



(5.19) 



By making use, first of Eq.(5.20) and then of Eq.(5.6), one can show that the primary con- 
straint surface V is preserved under the infinitesimal symmetry transformations, at least on-shell. 
Indeed,^^ 



[QT,Pm-fm{q,Pa,t)}p.B. 



5pr, 



dpa 



dq™- 



A. { 9l\ dL 1 \ 
dt [dqB J dqB \ j 



+ 5q^^ 



+ 



dq°- 



d / dL 
dt \dqB 



dL 
dqB 



Furthermore, we notice that the identity 



\dq-) dpa ■ 



dfm 
dpa 

(5.21) 
(5.22) 



^^The interplay between symmetries and conserved charges in the Lagrangian vs total Hamiltonian formalisms is 
briefly discussed in various exercises of the book [7] as particular cases of very general results on the elimination of 
auxiliary fields (mentioned in Footnote 22). More precisely, see e.g. Exercises 3.17, 3.25, 3.28 and 18.15 of [7]. An 
analogous derivation of such results for the extended Hamiltonian formalism should follow the general procedures 
introduced in [14]. As mentioned in the introduction, in the present text we prefer a more direct and pedestrian 
approach. 

^^In each line of (5.21) and (5.24), the velocity should be replaced by v{q,pa, g™, t) similarly to (3.27) and the 
expression is independent of the acceleration via cancelation due to (2.37), as it must be. 
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combined with the expression of momenta (3.4) when A = m, 

dL{q,q,t) 



5^" 



fmiq, q, t) , 



leads to 





fdL\ 


dL ' 


_dt 







(5.23) 



= 0. 
(5.24) 



5.1.3 Phase space variables 

As discussed in Sec.4.4, after solving the constraints, all the functions q'^ = (— q, t) 
are completely determined in terms of time and the phase space variables, together with the local 
free parameters v^{t) , as in Eq.(4.45). Therefore, all the velocitites have been removed so that the 
only independent variables are phase space variables. In other words, in this sense we are again 
working in the genuine Hamiltonian formalism. Substituting the general solution given in (4.45) 
into (5.16) reduces the infinitesimal transformations and the total Noether charge to be functions 
of {p,q,t): 

M^iP,q,t) := Sq^{q,pa,q"'ip,q,t),t) , 

Apa{p, q, t) := 5pA {q,Pa, q^{p, q, t),t) , (5.25) 

QT{p,q,t) := QT{q,Pa,q"'{p,q,t),t) = Aq'^pA- dK {q,pa,q"'{p,q,t),t) , 
which satisfy 

[Qt, /}p.B. = -A/ - Ag™}p.B. (Pm - Uq,Pa, t)) , 

[Qt,Pa}p.B. = -ApA - \pA,Aq''''}p,B.(^Pm " fm{q,Pa,t)^ ■ 



(5.26) 



The total Hamiltonian is equal to Ht :— v'^Pa — L{q, v, t) according to (4.47), where the 
expression of given by (4.45) is substituted. Now, we get on the primary constraint surface V 
that 

[Qt,H^}p.b. - [QT,v^}....PA-v^ApA+Aq^^^^^l^-[^^^^ (5.27) 
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due to (5.26). But the expressions of the momenta (3.4) show that the sum 



[ Qt. V [pa j ^[Qt:V j , (5.28) 



vanishes on the primary constraint surface V . Therefore, 
[Qt, Ht}p.b. 



d5K{q,v,t) _ dSq^{q,v,t) dL{q,v,t) \ A dL{q,v,t) 

QqA QqA Q^B j^'^l QqA (5.29) 

d6K{q,v, t) _ d6q^{q,v,t) dL{q,v,t) 
dt dt dv^ 

d{SK) d (Sq^) 



dt dt 



PA 



where we made use of Eq.(5.6) to get the third line and of Eq.(2.35) to obtain the fourth line. 
In terms of the very definition of the total Noether charge (5.25), we have thus shown that the 
Noether charge is conserved on the primary constraint surface, 

[Qt,Ht}p.b. + ^^0, (5.30) 

It is worth noting that this result is off-shell and parallel to the off-shell invariance of the action 
under the symmetry transformation. Of course, on-shell dQT/dt = [Qt, Ht}p.b. + — 
to be compared with the on-shell conservation (2.42) of the Noether charge. Another way of 
expressing (5.30) is to say that [ Qt, Ht}p.b. + is a linear combination of the constraints. In 
other words, the total Noether charge generates a transformation which preserves the Hamiltonian 
on the primary constraint surface. 

Furthermore, from the last formula in (5.18), not only Qt but also Qt preserves the primary 
constraints on-shell as in (5.24). More precisely, 

[QT,^m{q,Pa,t)}p,B.-^^ on-shcU. (5.31) 
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Therefore, the time evolution of this condition also vanishes on the primary constraint surface 
and on-shell, i.e. 



dt 



[QT,(pm{q,Pa,t)} 



P.B. 



^ d 

[Qt , <Pm}p,B. > ^t\p^ +-g^[QT, <Pm}p,B. - > On-shcll. 

(5.32) 

For the secondary constraints which essentially stem from [0^ , Ht}pb + ^If". we notice 



(5.33) 



P.B. 



and hence, from (5.30) and (5.32) we deduce that not only the primary constraints hut also the 
secondary constraints are preserved on-shell by Qt if [ Qt, Ht}pq + corresponds to a first 
class constraint. In this case, Qt is first class on-shell. As we see in the next subsection, the 
condition further implies that Qt preserves the solution space too. 



5.2 Symmetry in the total Hamiltonian system 

In this subsection, motivated by the results in the previous subsection where we studied the gen- 
eral properties of the total Noether charge which originates from a symmetry in a Lagrangian 
system, we discuss the symmetry in the total Hamiltonian system directly without referring to 
any Lagrangian system. In order to make the analysis concise, we introduce a single letter, x^^, 
1< M < 2J\f, to denote both p and q, 

x^ = q^, x^+^ = Pa ■ (5.34) 

We define a 2J\f x 2J\f constant non-degenerate graded skew-symmetric matrix by 

tMN _ _ _f_U#M*N T^M _ ( (-l)#^5^s 



which gives 



In particular, [x^,G]p^ = J^^SnG, where Sn = 



[F,G}B_B. = iLj^''^- (5.36) 
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5.2.1 Definition of symmetry transformations 

Now we define a symmetry of the total Hamiltonian system as a coordinate transformation on the 
jet space that 

1. depends on the phase space only, (i.e. it should not depend on , , etc), 

_^ x"'{x,t); (5.37) 

2. preserves the symplectic structure 

■>"" = ^ J"' %r = K. -'"}p.B. ; (5-38) 

3. takes any physical solution to another which means the preservation of both the on- shell 
relations and the constraints. Namely, if x{t) is a solution of the time evolution governed 
by a total Hamiltonian Ht{x, t; v, w) , given by (4.55), 

i:^ = [x^, Ht{x, t; V, w)}^^^ = J^'^'dNHrix, t; v, w) , (5.39) 

then so must be x'^'\x{t),t) for the same total Hamiltonian, up to the change of the local 
parameters and w^^' {t). 



X 



- x^'dNx"' + dtx"^ = j'^^'d'^HTix', t; v', w') . (5.40) 



Furthermore, such a symmetry must preserve the constraint surface V on-shell, 

(f)hix',t) c\{x,t)(f)g{x,t) on-shell. (5.41) 



Infinitesimally, the second requirement (5.38) reads 

Jm'MSx"^) = (-!)#-#- J^i9^(5x^) . (5.42) 

In other words, the 'super' one-form 5xm '■= Jm\^^^ closed, hence exact. Therefore, there 
exists a "generating function" Qh{x) on the phase space such that, c.f. (5.26), 

Sx^ =[x'',Qh}p^ . (5.43) 

Conversely, any such transformation leaves the symplectic structure invariant. 
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In comparison to the above definition of symmetry transformations, which refers to a specific 
given total Hamiltonian, a "canonical transformation," x x'{x,t) is defined as a coordinate 
transformation on phase space such that for every total Hamiltonian HT{x,t) there must exist 
another (not necessarily equal) total Hamiltonian H!p{x', t) obeying 

x'^ = x^'dNx'^ + dtx'^ = J^^'d'^H^rix', t) . (5.44) 

As the time evolution is generated by the total Hamiltonian, x^ = J^^dxHTix, t) , the condition 
(5.44) is equivalent to 

[x'^, x""}^^ &^Ht{x, t) + dtx"" = j'^'^&^HUx', t) , (5.45) 

since 

The analysis of Eq.(5.45) leads essentially to an integrability condition on the left hand side for 
arbitrary Ht{x, t), in order to be consistent with d'j^d'j^H'rp = [—1)*'^*^ d'jqd'j^Hlp. Namely with 
the notation x'j^ :— J]^\x'^, the integrability condition reads 

{cT^ [x^, x"<}^^)d'^HT + [x'^, x'^]p,sd'AHT + ^MX^'OKd^x'^ 

= {-lf^^*-{d'^[x'^,x'^]^^)d'^HT+[x'^^^^^ 

(5.47) 

This must hold for arbitrary Ht{x, t) and hence we have three independent relations: 

[x'^, x'^]^^) d'^Hr = a; [x'j^, x'^]^^ &kHt , (5.48) 

{[x'n,x'^}p.b) d'AHT = (-!)#-#- WM^x'^'j^^d'AHr, (5.49) 

ffj^x^'dKOtx'^ = {-l)*^*''d'^x''dKdtx'j^ . (5.50) 

Firstly, the second relation (5.49) with the quadratic choice Ht — x'^x'^ shows that [x'^, x'^}p^ 
is proportional to or 

K,a;'^}p.^. = /(a;,t)J^^. (5.51) 

Secondly, Eq.(5.48) with the linear choice Ht = x'^ further reveals that /(x, t) is independent 
of X , i.e. 

dMf{x,t)^0. 
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Finally, the last relation (5.50) shows that there exists a bosonic function $(a;',t) satisfying 
dtx'j^ = d']^j^{x', t). Using this and from (4.3), (5.51) we note that the explicit time derivative of 
f{t) vanishes as 



P.B. 



(5.52) 

= 0. 

Thus, from Eq.(5.51) one notices that canonical transformations leave the symplectic structure 
invariant up to a constant [x'^ , ^'^}pb ~ constant x J^^. Shortly, up to rescalings, canonical 
transformations are symplectic transformations. 

Finally we note that, if we require the preservation of the Z2-graduation and of the reality 
properties {5x^y = 5{x^'^), we may set Qh to be bosonic and Hermitian Qh — Q\j, so that 

{Sx'^y = [x^\ Qh}^^ = 5{x^^) . (5.53) 
5.2.2 Criteria for symmetry generators 

In order to clarify the criteria for the generating function Qh in (5.43) to meet the remaining con- 
ditions (5.40) and (5.41) as to be a symmetry generator, we investigate the infinitesimal version 
of them which are given by'^^ 

[ 5x^, Ht]p^ + dtiSx"") = J^^ (dx^^dLd^Hr + {d^(l>Y') 5v' + {dN<l>H)<l>h'Sw^'') 

"^■(5.54) 

and for 1 < /i < M + M', 

5x'^dM4>h = [x"^, Qi/}p.B. dM4>h = [Qh, 4>h}p.B. ~ on-shell. (5.55) 



"'"From the Leibniz rule of the Poisson bracket (4.4), it is worth to note an identity for an arbitrary function F{x, t). 



51 



The latter simply implies that Qh must be first-class on-shell. The former condition (5.54) must 
hold for arbitrary solutions of the total Hamiltonian. In particular, it should hold at the initial time, 
say to, at which the initial data xq can be taken arbitrarily. Thus the condition (5.54) should be 
interpreted off-shell, and the general solution 6x^{x, t) of the partial differential equation (5.54) 
may lead to all the symmetries in a given total Hamiltonian system. Rather, we translate the 
condition (5.54) in terms of the symmetry generator Qh, as done in (5.55), 



[x^, Qh} , Ht]^^ = - [Qh: K, Ht} + [x^, (l)l''6v' + '^Mh'Sw'^''' - dtQn} 



P.B. 

(5.56) 



Due to Jacobi identity this is equivalent to 

, [Qh, Ht}p,h. + dtQn - cl>]:''6v' - '^Mh'Sw^'''] = . 

J P.B. 

This condition should hold for every x^, I < M < 2J\f. Therefore, 

[Qh, Ht}p,s, + dtQn - (t>Y'5v' - iMw^w^""' = f{t) , 

where f{t) is an arbitrary time dependent, but phase-space independent function. This function 
can be removed by a redefinition of the generator. 



(5.57) 



(5.58) 



Qh 



Qh+ / dt'f{t'), 



(5.59) 



as the shift has no effect on the symmetry transformation, 6x = [x , Qh} p.p.- 



Thus, the necessary and sufficient condition for an on-shell first class quantity Qh{p, Q, t) to 
be a symmetry generator of a given total Hamiltonian Hxip, q, t; v, w) reads 

[Qh, Ht}ph. + ^ = + iMh'Sw^"' , (5.60) 

which is consistent with (5.30). The usual Hamiltonian version of the Noether theorem in uncon- 
strained systems states that any conserved charge (dQn/dt = 0) is a symmetry generator. The 
formula (5.60) is the corresponding generalization to constrained systems. 



5.3 Solutions 

• Every quantity, which is first-class and conserved on-shell, is a symmetry generator 
Indeed, the fact that Qh remains constant under the time evolution reads 

^^[Qh,Ht}ph. + ^ = 0, (5.61) 
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which is stronger than (5.60). 

• The total Noether charge Qt (5.25) originating from a Lagrangian system is a symmetry 
generator if [Qt, Ht}p^ + corresponds to a first class constraint. In this case, the 
full expressions of the right hand sides of (5.25) should correspond to Sx'^ and it follows 
that Qt is first class on-shell, as we saw in Sec.5.1. 



Static examples: 

In the static case, there is no explicit time dependence in the constraints (f)h{p, q) and the 
quantity H'{p, q) defined in (4.56). The time derivative of any first-class constraint (t)Y^ is 
then first-class too, as shown in (4.69). 

- In the static case, the total Hamiltonian itself corresponds to a symmetry generator 

at at 



- If there is no secondary first-class constraint in the given system (so that all the first 
class constraints are linear in (f)]^^), then 

= [(l>l^\ Ht}^^ = <t>fT\ . (5.63) 

Hence <^}^^s\t) corresponds to a gauge symmetry generator with arbitrary time de- 
pendent functions s^{t), i.e. 

Qh = (t>Y's\t) , (5.64) 

satisfying (5.60). 

- Alternatively if the time derivative ofcj)]^^ is quadratic in the constraints 4>h (and hence 
first-class) 

^l^' = Ht]p^ = IMh'T^''^ , (5.65) 
then again (f)l^^6''{t) corresponds to a gauge symmetry generator, 

Qh = <l>l''e'{t) . (5.66) 
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- Combining the above two cases we have more general solutions. Namely, if the time 
derivative ofcfP"^^ is a sum of terms linear in and quadratic in (j)^, 

4>':' = Ht]^^ = (t>l''T\ + |0,0,.r'''^'„ , (5.67) 

then (pl^^e'^it) corresponds to a gauge symmetry generator, 

Qh = (Pa'e'^it) . (5.68) 

For example, we consider the Lagrangian L{x, y, x, y) = \ e^x^, whose equations of 
motion leave y arbitrary (so y is pure gauge) but fix x in time x = xq. This Lagrangian 
produces the Hamiltonian H = \ e~^p^, one primary first-class constraint (j)^^^ = Py, 
one secondary first-class constraint px, and the total Hamiltonian Ht — \ e^-'pl + 
Pyv{t) . This is a counterexample to Dirac's conjecture (see e.g. the subsection 1.2.2 
of [7]) because the secondary first-class constraint p^ does not generate any gauge 
symmetry as x is fixed by the equations of motion. However, the time derivative of 
the primary first-class constraint satisfies Eq.(5.65), 

Py = \py, Hrjps. = I (^'^Pxf , (5.69) 
and generates arbitrary shifts of the pure gauge variable y . 



- A linear combination of the primary and secondary first-class constraints, 4>]^^ and 
(pl^^, can be a gauge symmetry generator, 

QH^cj>l^'e''{t) + cl>l^'e%t), (5.70) 

if the local functions, s^{t) ande^{t) satisfy, with (4.70), 



dt 



+ T%e\t)+T\e\t) = Q. (5.71) 



For example we consider the Maxwell theory of which the Lagrangian and the Hamil- 
tonian read"^^ 

C — — i F Fi^^ 

H^H' = j d^-'x [\FijFij + \WW - Aadilli\ , ^5 ^72) 

Ht^ j d^-'x [\FijFij + \WW - A^diW + v{x)Il''] . 



' This example is also handled in the section 1 9 . 1 . 1 of the reference [7] . 
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The gauge symmetry is one example of (5.70) and (5.71) as 



I 



dx^-^ [e{x) d0 - doe{x)Ii''] , 



(5.73) 



where 11^ = are the gauge invariant canonical momenta for A^, and 11° is the 
primary first-class constraint, while dilT is the secondary first-class constraint. There 
appears no other constraint. 

- In the extended Hamiltonian formalism of Sec.4.6, every first-class constraint corre- 
sponds to a gauge symmetry, if the system is static. Namely (j)}^^e'^{t) + 0jfV*'(t) is a 
gauge symmetry generator with arbitrary time dependent functions, e^{t) and (t). 



As one can see, the generators of local (i.e. gauge) symmetries are linear combinations 
of the constraints, therefore they vanish weakly in contradistinction with the generators 
of global symmetries. In this sense, only global symmetries lead to non-trivial conserved 
charges. 

5.4 Dynamics with the arbitrariness - gauge symmetry 

We remind the reader that the key philosophy we stick to is that different choices of the local 
gauge parameters correspond to different total Hamiltonian systems, which nevertheless should 
be taken equivalent, describing the same physics (see Subsection 4.5.2). 

A somewhat less drastic - though equivalent - perspective is to consider only one total Hamil- 
tonian throughout the time evolution, with a single set of local gauge parameters. The local 
functions should be continuous all the time but infinitely differentiable, i.e. C°°, only piecewise 
in time. This discontinuity in the derivatives corresponds to changes of local gauge parameters 
at different times. As long as the local functions v'^{t) are continuous, one can change them arbi- 
trarily at any moment. The continuity guarantees the continuity of the first order time derivative 
of the dynamical variables (p, q). However, in the Hamiltonian dynamics, there is no reason to 
require the continuities for the higher order time derivatives. 

Explicitly, expressing the dynamical variable, F, at a future time, t + dt, as a power expansion 



it). 



(5.74) 
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of dt around the present time, t, we have 



F(t + dt) ^ F{p{t'),q{t'),t%,__ 



4+dt 



^F + dtF+ldt^F + 0{dt^) 



r. , fdF 

d^F 



F,H, 



+ \de 



F + dt 



F,Ht] 



+ 2 



dF 
'dt' 



+ 



P.B. 



P.B. 



+ Idt^ 



' dt 

[F,Ht},Ht} +0{dt'), 

J P.B. 



[F, Ht} , Ht 



P.B. 



where we have assumed that F exists or v'^{t) is differentiable, and we have set 

F^F^dt—^\de^^ 



(5.75) 



+ 1 dt^ ^H' + + \ dt 



dH' d(Pj'' , ^^^^dv' 

\ -^v + 6 — 

dt dt dt 



(5.76) 



The coefficients v^{t) are completely arbitrary and at our disposal. We recall that different 
choices of the local parameters mean different total Hamiltonian systems, which nevertheless 
should be regarded equivalent, i.e. describing the same physics. For two different choices of the 
coefficients, say v and v + Av, the dynamical variable at the future time differs by 



AF{t + dt) 



dt 



' ) P.B. 



^dt' 



+\dt^ 



+\dt^ 



F, 



501 



1st 



dt 



P.B. 



Av' + Idt" 



A dAv' 
ip.B. dt 



F, 01^*} , i/'}^^ A^;^ + ^dt^ [[f, //'1 , 0''*}^^ ^v' 
F, (j)]'' I , (PY' I (v'Av^ + v^Av' + Av^Av') 

) ) P.B. 



(5.77) 



+ O {dt\ Av) . 
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Thus, the leading order in the difference appears at the first order in dt or the velocity, when 
Av{t) 7^ 0. Namely, different velocities for the same initial configuration can still correspond 
to the same physical state. This may"*^ correspond to the choice of the Noether charge, Qt — 

(/)^"*£»(i) such that e''{t) = and i'(t) = Av'' at time ' t '. 

On the other hand, if at time t one has Av{t) — 0, then the time derivative (p, q) of all the 
dynamical variables are the same in the two different total Hamiltonian systems, and the first 
nontrivial difference appears at the order of dt"^ or the 'acceleration'. 



Again, this may correspond to the choice of the Noether charge, Qt — 0j®*£'(i) such that e^{t) — 
0, e'{t) = and e'{t) = ^ at time t. 



'^^For Qh = 0i*'*e*(t) to be actually a symmetry generator, meaning it preserves the solution space, {p{t), q{t)}, 
some extra conditions should be satisfied as (5.67) or (5.70). 




(5.78) 
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6 Dirac quantization for second class constraints 
6.1 Dirac bracket 

On the 2A/'-dimensional phase space with variables and (1 < A, -B < A/"), we consider a 
set of functions q) (where 1 < s < dim{ps} < 2jV) such that the following supermatrix is 
non-degenerate, 

Qst [Ps,Pt}p.B. = -i-lf'*'nts , (6.1) 

or its inverse exists, 

= ^ {n-'y'ntu = 5\. (6.2) 

We note 

(n-Y = -(-!)#«#*+#«+#* (n-y^ , (6.3) 

and for an arbitrary quantity, F, 

We define the "Dirac bracket" associated with the functions ps as, 

[F, GU,^^ := [F, G}p.B. - [F, Ps}p.b.^-'''[pu G}p.b. • (6.5) 
Some crucial identities follow. We first note that for an arbitrary object, F, 

[ps, i^joirac = , [F, pjnirac = . (6.6) 

This property is the raison d'etre of the Dirac bracket. It means that one may impose ps = 
either before or after computing the Dirac bracket, whichever one prefers. Just like the Poisson 
bracket, the Dirac bracket satisfies the symmetric property, 

[F, G},,^ = -{-l)Mo [G, F},,^^ , (6.7) 

and the Leibniz rule, 

[F, GX}o.ac = [F, G},,^^ K + {-l)*-*oG [F, K}^._ , 

(6.8) 

[FG, K}^_ = F [G, K}^,^^ + (-l)#o#- [F, K}^_ G . 

Moreover, from the Jacobi identity for the Poisson bracket, (4.6) and (6.3, 6.4), one can verify the 
Jacobi identity for the Dirac bracket after some tedious calculations, 

m G},,., , i^},,., = [F, [G, Hj^^^J^^^^ - {-l)Mo [G, [F, H}^._}^^^^ . (6.9) 

In mathematical terms, one says that the Dirac bracket obeys to the axioms of a graded Poisson 
bracket. 
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6.2 Total Hamiltonian dynamics with Dirac bracket 

One can prove by contradiction that the Poisson bracket between all the second class constraints 
is non-degenerate. For constrained systems, the Dirac bracket is defined for all the second-class 
constraints by identifying ps with (j)^^'^ . It is convenient to let the Dirac bracket governs the 
dynamics rather than the Poisson bracket, 

l£(^=lf.H,}„, +f . (6.10) 

because in such case one may impose the second-class constraints before computing the evolution 

of the system. It follows that one may omit the second-class primary constraints while adding to 
the Hamiltonian or to the Noether charge in the definition of the total Hamiltonian (3.30) or the 
total Noether charge (5. 17). We note that 

[F. Ht}p,^, - [F, Ht},._ = [F, <t>':^}p.Bn-'^\<t>r\ Ht}p.b. ■ (6.11) 

Hence, as long as the second-class constraints have no explicit time dependence, the Poisson 
brackets [0^"*^, Ht}p.b. in the right-hand-side vanish on the constraint hypersurface V . In such 
case, the dynamics with the Dirac bracket and the other with the Poisson bracket, are identical on 
V. Namely both reduce to the same Lagrangian dynamics. 

• First order kinetic terms 

In most of the cases, the momenta Pa for the fermions are linear in the spinor field V'°> 
resulting in the primary second-class constraints, 

(l>^a^ ■■= Pa - L^pll^^ , L^p = Lpa , [02-d, (^f^}pB. = 2L^p . (6.12) 

Using Eq.(6.5), the Dirac bracket of an unconstrained bosonic system coupled with fermions 
read 

bosons fermions 



where the factor 1/2 comes from the factor two in the last equation of (6. 12). Notice that it 



is important in Eq.(6. 13) to treat the partial derivatives {q^, pb, t/j", ^l'^'^) as the independent 



variables rather than pi,, p^), and hence 

= , [r, V^n_ = -1^-^"^ , [Pa, V^^}_ = , etc. 

(6.14) 
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The last expression should be compared with [pa,'^l^^}pg = Sa^- This 'halfness' of 
the Dirac bracket is not related to the fermionic character of instead it is typical for the 
system with first order Lagrangians (such as the variational principle of the Hamiltonian 
formulation itself). If tjj is complex, then the 'halfness' of the quantization is 'doubled' and 
one recovers the standard naive rule of the canonical quantization for a Dirac spinor. 

• Quantization 

The quantization can be straightforwardly performed by replacing the Dirac bracket by the 
super-commutator with a factor —i, 

[ , lo^ac =^ , }, (6.15) 

which gives the standard convention, [q^tPb} — + i ^^b- The point is that, from 

the second-class constraints are central, even after the quantization. Therefore, one can 

simultaneously impose the second-class constraints on the physical states. This would not 
be possible if one had naively performed the correspondence rule in terms of the Poisson 
bracket itself. The second-class constraints should be represented by identically vanish- 
ing operators on the physical Hilbert space. In practice, this may be realized by solving 
explicitly the constraints in terms of some set of independent variables 

0r'^(a:^t) = O ^ x^ = f''{f) (6.16) 

and try to represent the algebra [y", y^'jDirac — g^^iv) on the Hilbert space of functions of 
the y's only. 

Although this way of quantizing second-class constrained systems looks pretty straight- 
forward and conceptually clear (one imposes all the constraints), in most practical cases, 
second-class constraints are most tedious because in general either we are not able to invert 
the matrix Qgt and the Dirac bracket is not known explicitly, so that nothing can be done 
at all, or we are not able to find a faithful representation of the Dirac bracket algebra."^^ 
Somehow surprisingly, first-class constraints are preferable because there is an algortihmic 
- though involved and subtle - way to quantize the theory in terms of the Poisson bracket 
(which is easy to represent). 



"*^More comments on the quantization of second-class constraints can be found in the section 13.1 of [7J. Notice 
that some systems admit only the Dirac method of quantization and not the so-called "reduced phase space" method 
[18]. 
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7 BRST quantization for first class constraints 

The BRST procedure is motivated through the Faddeev-Popov construction. Here we review the 
essential features of them in a self-contained manner. 

7.1 Integration over a Lie group - Haar measure 

For a Lie group G of dimension Ng , we parameterize its elements g by the coordinates 6*" (1 < 
a < Ng) of the corresponding Lie algebra of a basis {Ta}, 

g{e) ^ e''"^^ , [T,,T,]^iC:,n. (7.1) 

We also define a set of Ng functions C"(^i; ^2) from the multiplication, 

^(^1)^(^2) =^(C(^i, ^2)) . (7.2) 

From the Baker-Campbell-Hausdorff formula. 

In (e'^e^) — x + y + ^[x,y] + higher order commutators , (7.3) 

we obtain explicitly, 

(^(61, 62) = 9l + e'i-\ e^^e^C^l + higher order terms . (7.4) 
The left invariant measure for the integration over the Lie group G is denoted by 

Ng 

Vlq := vewLie) = H WLie) . (7.5) 

a=l 

By definition, it must satisfy the property of left invariance, i.e. for an arbitrary fixed element 
go ^ G and any function on the group F{g), 

j VLgF{g) = j VLgF{gog) . (7.6) 



Hence, the following identity must hold for any 6q and 9, 



w^{e) = det ( ^^^i^ ) (C(^o, e)) . {1.1) 
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Some simple choices like ^ = or = — ^ give explicitly 



WL{0)^WLiO) det 



89^ 



Wl{Q) det" 



1 /ac(M) 



(7.8) 



i?=0 



Similarly one can define the right invariant measure V^g — V9Wr{9) to obtain 



Wr{9) = Wr{Q) det 



dC{9,9^) 
d9^ 



Wr{Q) det" 



(7.9) 



Now we are going to show that both measures can be set equal. The chain rule for 9' :— 
C{9,v) gives 



dC{C{9,i9),9o) 



dC{9',9o 
89"^ 



dC\9,^) 



i?=0 



because 9' — 9 when v = 0. The associativity property explicitly reads 

C(^i, 92, 9s) - C(C(^i, 92), 9s) = C(^i, ({92, 93)) . 
Therefore, evaluating Eq.(7.10) at = — ^ > we obtain 



det 



dC{9,i9,-9) 



= det 



i?=0 



dC{9,9o 
89^ 



det 



dC\9,d) 



(7.10) 



(7.11) 



(7.12) 



From 



g{a9,^, -6)) = g{9)g{^)g{9)-' = e'^^^^o^ 



it follows that 
Consequently, 



Tea := g{9)Tag{9)-' =: {Me)a% , (7.13) 
C{9,^,-9) = i)\Mo)b\ (7.14) 



tT{Tan) = {Me)a\Me)b''iT{TcTd) . 



(7.15) 



Thus, as long as ti(TaTij) is invertible as a Ng x Ng matrix, {e.g. when Ta are in the adjoint 
representation of a semisimple'^'* Lie group) we have (detM^)^ = 1. Furthermore, from the 
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If it is compact, then one can further take tr(TaT(,) cx 5ab- 
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continuity at ^ = we disregard the possibility of det being minus one. Thus, we obtain from 
(7.14) 

«(^?^^)=detM. = .l. (7.16, 

Inserting this relation on the left-hand- side of (7.12) and setting Wl{0) = ^^^(0), one has 
shown that the left and right invariant measures may be chosen identical for the Lie groups where 

det (tr(r„n)) 7^0: 

Wnie) := Wd0) = Wnie) . (7.17) 
The corresponding measure is known as "Haar measure" [19] and satisfies 

J VgF{g) = j VgF{go9) = J VgF{ggo) = J VeWH{e)F{g{e)) . (7.18) 



In the case of the local gauge symmetry in field theories with the gauge group Q, the parame- 
ters 9"'{x) are in fact arbitrary local functions. We assume that there exists a countable complete 
set in the commutative algebra of local functions, 

{fn{.x)}, fn{x)f^{x) := d^^fi{x) , {fn{x)\fm{x)) = S^m , (7.19) 

where d^^ are the structure constants of the algebra with the pointwise product. Then we can 
write e^lx) = ^""/„(x) so that 

(7.20) 

As the structure constants for the set {Tan{x)} are independent of the coordinates x, the above 
relations induce a novel group ^ = | 90"'^) |> which is defined by the representation g{9{x)), 

with the parameters Namely, though the representation is given for some fixed coordinates 
system x, there exists an abstract group independent of the coordinate choice. We have 

9{k)m) = 9 {((01,02)) , C"(^\,^2) := {fn{x)\C{er{x),e2{x))) . (7.21) 

The group G is infinite-dimensional since it is a local (i.e. position dependent) group. 



63 



Now we are ready to straightforwardly apply the left/right invariant measure to this local 
group: 



j VgF{g) = J 



VgF{g) = / VeWH{e)F{g{e)) = / VgF{gog) = / VgF{ggo), 



Wh{0) = det 



.1 ( d C""{d,e ) 

Q^bn 



det 



i9=0 



Q^bn 



(7.22) 

The situation in gauge field theories is that F{g) is given by a spacetime integral of a functional 
F{g, d^j,g) which depends on the local group element g and its spacetime derivatives d^j,g, 

F{g)^ j d''xFig,d^g). (7.23) 

In the remaining of the paper, for short notation we drop the hat symbol and simply denote e.g. 

J VgF[g] = J VgF{g). (7.24) 

7.2 Faddeev-Popov method 

We consider a dynamical system where a finite-dimensional Lie group Q, acts on the dynamical 
variables which we denote collectively by For each element, g & define a map (i.e. here 
a gauge transformation), 

9 : ^^^g- (V.25) 

For later purpose, we write the successive gauge transformation in the following order, $ — > 
— > $3251- other words, one has a left action of Q on the space of dynamical variables. 

We introduce a Lie algebra valued functional /?.[$] = /i"[$]Ta of $ (which may depend on its 
derivatives as well). We assume that is non-degenerate under the gauge transformations as"^^ 



det 



(7.26) 



6»=0 



in which case it is called the "gauge-fixing functional". 



''^It is not necessary to require the non-degeneracy for 6^0. 
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For an arbitrary function B{h) of h"-, the Faddeev-Popov functional Bp p [$] of $ reads 



Sp.p.[$] := B{h[^]) det 



(7.27) 



By making use of the definition (7.22) of the Haar measure, one finds the following crucial iden- 
tity. 



jvgB^„{^g\ ^ jvgB{h[%]) det 



I 
I 
I 



Ve det 



= / ve det 



89'' 
89'' 



8-0^ 



V9 det ( ) B{h[%,e)]) 



= / VhB{h), 



where 



Ng 

Vh := JJd/i" 



(7.28) 



(7.29) 



a=l 



Thus, the integral of Bp p, over Q is just a c-number, being independent of $ or the choice of 
the gauge-fixing functionals provided that the integral domain for h is $ independent. 



Now we consider a gauge invariant functional, 
and further assume that the path integral measure is gauge invariant. 



det 



1, 



(7.30) 



(7.31) 



which is always satisfied for any dynamical system where the fields are in the adjoint representa- 
tion or there are equal number of fundamental and anti-fundamental fields, transformed by g and 
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g ^ respectively. 

The Faddeev-Popov path integral method prescnbes to multiply the gauge invanant functional 
integrand by the Faddeev-Popov functional. 



=0 

(7.32) 



J J^{^) — ^ Jv^ Bp.p.[$] = J'^^ J^{^)B{hm) det 
to satisfy 

JvgJ P$.F($)Bp.p.[$] - j'^aj V^J^{^)B^,,,[^g] = Q VhB{h)^ j V^T{^), 

(7.33) 

where we made use of (7.28), (7.30) and (7.31). Equivalently we have 

I P$J'^($)i3p.p.[$] = ( / VhB{h)] oc / V^T{^). (7.34) 



In gauge field theories, the volume integral of the gauge group is often divergent, and the 
Faddeev-Popov method [16] provides a regularization scheme for that by adding the Faddeev- 
Popov factor to the original gauge invariant action. In the path integral, the gauge invariant 
functional is given by the exponential of the gauge symmetric action, 

= e+^'5[*l , cS[$] = S[^g] . (7.35) 

Furthermore, in this context, the arbitrary function B{h) is reformulated as a Fourier transforma- 
tion, 

B = j DA;e+'(-^('=) + *^<''^"[*l) , (7.36) 

where the fields A;„ are said to be "auxiliary" and V{k) is proportional to the logarithm of the 
Fourier transform B{k) of the function B{h) (anyway, the coice of V is as arbitrary as the one of 
B). Moreover, the Faddeev-Popov determinant can be written by introducing a pair of fermionic 
scalar fields, called "ghosts", uoa and cj* (1 < a, 6 < Ng), 

det = / -DujVuj e+^'^<'^"^[*l-' , (7.37) 
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where we set"*^ 



(7.38) 

6»=0 



Note that uja and uj^ are not necessarily complex conjugate to each other. We further assign the 
ghost number, +1 for uj"-, —1 for 0;^, and for the other fields $ and 

Combining them all, the Faddeev-Popov action reads 

5p.p. [$, a;, a;, A;] = - V{k) + A;„r [$] + a;„A%[$]a;'' . (7.39) 



For example, in the Yang-Mills theory (on a curved spacetime), if we take the Lorentz gauge, 

h[A] := V^A^ (7.40) 

then we can set 

5p.p.[$,u;,(D,A;] = Sy^l^ + j dx^^/gtr(^ - V{k) + kh[A\ - V^ujD^uj^ . (7.41) 

Another example is the gauged Hermitian one-matrix model. Diagonal gauge choice leads to a 
Faddeev-Popov determinant which is nothing but the Vandermonde determinant [17]. 

7.3 BRST symmetry 

Remarkably, even after choosing a gauge"^^, the path integral still does have a symmetry related 
to the gauge invariance. Indeed, the Faddeev-Popov action (7.39) possesses a fermionic nilpotent 
rigid symmetry known as "BRST symmetry", after its discoverers, Becchi-Rouet-Stora [20] and, 
independently, Tyutin [21]. 

To discuss the BRST symmetry it is useful to note that the infinitesimal transformation asso- 
ciated with the Lie algebra element Sg — i ■d"'Ta is given by 



(7.42) 

s=0 



^^All these relations are valid for theories in the Minkowskian spacetime. For the Euclidean theories we only need 
to replace the factor, +i by —1, and consider the Laplace transformation rather than the Fourier transformations. 
"''For example, the choice V(k) =0 leads to a delta function to fix the gauge as = 0. 
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and, from (7.3), the Faddeev-Popov matrix A"b[^] transforms as. 



while transforms as 



s=0 



ds V de^ 



e=o, s=o 



C=o 



ds 



e=o, s=o 



(7.43) 



(7.44) 



With a fermionic rigid (i.e. independent of x) scalar parameter e, the BRST transformation 
reads as a 'gauge transformation' generated by i?" = eo;" or, equivalently by 

Sg^ieu^Ta, (7.45) 

so that 



5$ = — $ 

ds 



(/(sew) 



= , 

s=0 



(7.46) 



5a;" ^ -\s C^uj uj" , 5uja ^ -eka, 



where we also defined the transformation of the ghosts u and uj . 

In Yang-Mills theories the dynamical variables {$} = {A^, cf), ^, tjj} consist of a vector field 
A^, matter fields ijj, ip in the fundamental, anti-fundamental representations and matter fields in 
the adjoint representation,"^^ such that the standard gauge transformations are 

Ag^, = gAi^g-^ - idfj^gg-^ , 0g = g(pg~^ , i/jg = gi> , i/jg = i/jg'^ , etc. 

(7.47) 

Explicit expressions for the infinitesimal BRST transformations read [22] 

SA^ = eD^oj , (50 = i[ea; , 0] , 5ip = ieuip , Sip = —itpeuj , (7.48) 



our analysis, the matter fields (p, ip and tj) can be either bosonic or fermionic. 
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where we set 

u-.^u'^Ta. (7.49) 

We introduce the BRST charge Qbrst which is the fermionic generator of the BRST transfor- 
mations satisfying 5^ — £[(5brst , More explicitly, 

[Qbrst , ^/i] = +D^U , [Qbrst , (f)} ^ +i[uj , (j)} , 

[Qbrst , -0} = +«<^'0 , [Qbrst , -0} = Ti'ipi^ , 



(7.50) 



{Qbrst , '^"} — ~2^bc'^ ) {Qbrst , '^a} — ~ka , 



[QbRST ) , 

where the sign, =p depends on whether ip is bosonic or fermionic. 

In particular, the transformation was designed to satisfy {Qbrst , A^^f^Jo;''} = in order for 
the Faddeev-Popov action (7.39) to be BRST closed, and such that 

{Qbrst, Co"} = +ia;^. (7.51) 

The BRST charge Q BRST increases the ghost number by +1. 

Because — 0, the finite transformations are then given by 

e-dQBRST = l + £adQ3^3,, (7-52) 
and, with g :— e*^'^ = 1 + ieu , each field transforms to 

$g , Ug^ gu^u^ ieiJ' , Uga = cDa - ska , kga = ka . (7.53) 
Further under the successive transformations, e'^^'^'^'^BRSTe^^adQgj^g^ ^ each field transforms as 

^ ^92 ^9l°92 93 1 

uj — > ujgr^= uj + ie2Uj'^ — > ujg^+ ie2{ujg^)'^ = uj + i{ei + e2)uj'^ = ujgs , 

(7.54) 

<^a ^ ^92a — ^a — ^ika <^gia " ^ikg^a = — (^1 + ^2)^0 = ^93a , 

"'a ' "'92a — "'a ' "'gia — "'a — "'gsa i 
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where we set 

5(2 := 1 + ie20Jg^ = 1 + ie20J + £1^2'^^ , ■= e'^^^^^^^'^ = 1 + i{ei + e-2)uj = gio g'2- 

(7.55) 

Thus, we note 

g£iadQgj^g^g£2adQgj^g^ ^ g(£i+£2)adQgj^g^ ^ (7.56) 

and hence, BRST charge is nilpotent, 

(ad^BRST)^ = , or [Qbrst , [Qbrst , Field}} = . (7.57) 
The nilpotent property can be directly checked from (7.50), using the Jacobi identity. 

From the gauge invariance of the original action it follows that is QaRST-closed, 

gBRST,5[$]] =0. (7.58) 

Hence, the Faddeev-Popov action reads as a sum of QBRST-closed and QeRST-exact terms, 

5p.p.[$,a;,a;, k] = + {Qbrst , i^aV'ik) - uOah'^m] , (7.59) 

where, without loss of generality we have shifted the arbitrary function of the auxiliary field, 
V{k), by a constant in order to satisfy V{0) — and to write it as 

v{k) := kaV'ik) = - jgBRST , ojaVik)] . am 



In particular, for the Yang-Mills theory with the Lorentz gauge (7.40), one can rewrite the 
whole action (7.41) as 

SrMm + j d2:-° v^tr( - V{k) + kV^A^' - V^uD^u) 

(7.61) 
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7.4 Hodge charge and a number operator 

The Hodge charge is defined to be fermionic and acts only on the auxiliary fields ka as 

[Qnodge , ka] = -Oa , [^Hodge , others} = . (7.62) 

Hence, it is nilpotent, 

(adQHoage)' = 0, (7.63) 

and satisfies 

{Qnodge , Qbrst} = N^,k , (7.64) 

where N^^^k is the number operator counting the total number of uja and kf, fields, 

[Nu>,k , (^a] = , [Nu,,k , ka] = ka , [N^,k , otheis ] = . (7.65) 

Both of Qbrst and Qnodge do not change the total number of Oa and kb fields since they do not 
annihilate the quantity such as {Qnodge, '^o} = . One can show straightforwardly that 

[QbRST , Ar^,fc] = , [QHodge , iVc^,fc] = . (7.66) 



Consider a QeRST-closed quantity T, 

[Qbrst, T} = 0. (7.67) 

The latter condition is extension of the gauge invariance condition for functional depending on 
$ only. We can decompose it as a sum of N^^^k eigenstates, 

oo 

T = ^ Tjv , [N^,k , Tj,] = NTn . (7.68) 

iV=0 

From [Qbrst , Na!,k\ — one can deduce that 

[Qbrst, T^} = 0. (7.69) 
Therefore, any QeRST-closed quantity can be written as 

<9brst,t} = ^ T[$,u;,cI;,A;] = To[u;,$]+ [Qbrst, T[$,u;,u;,/c]} , (7.70) 
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where 

(7.71) 

In other words, the cohomology group of the BRST charge is trivial at non-vanishing grading 
Nuj,k- This reasoning is a particular example of a standard procedure for computing cohomology 
groups. In mathematical terms, the Hodge differential Qnodgc is called a "contracting homotopy" 
for the number operator N^^^k with respect to the BRST differential Qbrst • More materials on 
general cohomology groups can be found in the introduction for physicists to graded differential 
algebras in the chapter 8 of [7]. 

The importance of the BRST cohomology sits in the general theorem that the BRST coho- 
mology group at ghost number zero is isomorphic'*'^ to the algebra of observables of the theory. 
In other words, for the quantum theory it is isomorphic to the physical spectrum. Apart from this, 
as was roughly shown in the above example, the importance of the BRST formalism in gauge 
theories is that it allows to write a gauge-fixed path integral and to make sure that the final results 
are independent of the choice of gauge (because the corresponding terms are BRST trivial), see 
e.g. [23]. Nevertheless, the BRST cohomology is of high interest already at the classical level, 
as explained in the report [11] where some examples of applications are given. For further dis- 
cussion on the 'Hamiltonian' BRST formalism presented here, we refer to the chapters 9 till 12 
and the chapter 14 of [7]. The 'Lagrangian' BRST formalism of Batalin and Vilkovisky [24] has 
the advantage of being entirely general (it includes the case of gauge algebras that close only 
on-shell) and covariant (since it is Lagrangian). The Batalin- Vilkovisky (also called "antifield") 
formalism is explained in the specific reviews [25] and in the chapters 17 & 18 of [7]. 



Qbrst, To[u;,$]| = 0, 



'The reader may consult the section 11 . 1 of [7] for a proof and for more conmients. 
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Appendix 



A Some proofs 

Here we present the proofs of some facts discussed in the main body of the text. 
• Eq.(2.12). 

Proposition. An arbitrary function on the jet space t) is a total derivative if and only 
if its Euler-Lagrange equations vanish identically, 

^Z(,„,t) = o ^ FM = '-^. (A.1) 



Proof: The proof of the necessity '<^' is straightforward from (2.10). In order to show 
the sufficiency, we filter the set of all functions on the jet space by sets J^n, 

:= {F{qt, t) , < n < iV} , TV = 0, 1, 2, 3, • • • . (A.2) 

We prove the sufficiency by mathematical induction on N : 

- When iV = 0, the left hand side of the claim (A.l) implies that the function 
depends, at most, only on the explicit time, t, being independent of q^, i.e. F{i). 
We can simply set K{t) = /J df F{t'). 

- Now we assume that the converse is true for any < < M, and consider the 

/d V 

case, N = M. It is useful to note that if F e J^n, then — 1 F e J-'n+v. and 



its only dependence on qn+u appears as 



dtj 



/ d \" dF 

^ ^ lN+n 'q^ + ^N+n-l , C^AT+n-l G ■ (A.3) 



Hence for F e J^m, from 

M 



we first note that F e J^m is at most linear in qM, i.e. 

F{qn: t) = qtiFA{qn, t) + Om-1 , Fa G Tm-i ■ (A.5) 
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Consequently, 



= 



5q^ 

M / 1 \ M-1 



(A.6) 

Thus, 

^-(-l)--lp^-O, (A.7) 

SO, by the usual Poincare lemma in the space of one-forms -?^a(?m-i)' there exists 
a function K{qn, t), such that 



i^A = ^^S^ , i^(gn, t) e . (A.8) 



dK 

Finally, if we define F := F — , then 

at 

5F' 

^ = 0, F'e^M-i. (A.9) 

Thus, from the induction hypothesis, F' is a total derivative, and hence so is F 
itself. 

This completes our proof. □ 
Eq.(2.16). 

Using Eq.(2.11), direct manipulation shows the following chain of identities, for an ar- 
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bitrary function F{qn,t) on the jet space and for m > 1 , 



n=0 



dtj \dq^ 



• Eq.(2.27). 



n=0 



= E 

n=0 
oo 

= E 

n=0 

oo 

= E 

n=0 



E(-^ 



dt V 9?;^ 



+ 



dtj 



■-T 

dtj 



d 



dt V dq^ 

dq^- 

dq 



m—l 



idF 



dq^ dt 



4- 



(A. 10) 



We show the relation (2.27) by induction on the power of s. We assume that the following 
relation is true up to the power A?^ — 1 in s. 



d 



(A. 11) 



which is clearly true for = 1, as g;^ = if s = 0. Now differentiating the left hand side 
with respect to s, we get, up to the power iV in s, 

d_ ( ,.n. d 

ds 



i-Affi(im,t)-g^q^ 



= fFiqm.t)-^ [f^{qk.t)^t 



d 



fl{qm:t)^mk:t) + 0{s") 



fri<lp^tmi^fniqm,t)+0{s^) 



^ dq^ d 



dqf dq[ 



ds dq^ 
d 



^Jn{qm,t) + 0{S 



N\ 



(A. 12) 
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Thus, the relation holds up to power N, and this completes our proof. 
Eq.(4.67). 

Proposition. For an arbitrary quantity, F{p,q,t), the two actions, namely taking the time 
derivative and taking the restriction on V, commute each order. 



Proof: Using the coordinates, {x, 

F{x,(f),t) = F{p,q,t). Then 



for the whole phase space, (4.42), we first define 



d 



;dF . OF OF 



^^dF{x,0,t) ^ dF{x,0,t) 



dx^ 



dt 



(A.13) 



□ 



-F(a;,0,i) 



Eq.(4.66). 

Proposition. The Poisson bracket is independent of time, 

[ ' }p.B. ~ [ ' }p.s. Ito 

or 



(-1) 



d d 



d d 



QqA 



(-1) 



d d 



d d 



(A. 14) 



(A.15) 



Roughly speaking, the Poisson bracket is independent of time, i.e. preserved, because time 
evolution is a symplectic transformation generated by the Hamiltonian. 
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Proof: We show the proposition for a set {f *(t)} of local functions which are piecewise 
infinitely differentiable. Then the time independence holds globally, since {p, q) are 
globally continuous. A direct manipulation gives 



[ , } 



P.B.I to 



d d 



d d 



d 



dq"^ 



d 



d 



P.BAto dqB QqA 



d 



P.BMo dpB 



(A.16) 

Obviously, the equality we want to show holds for the zeroth order in t — to- Now we 
suppose that it holds up to the order [t — toY, so that up to the power k. 



W^.Pb] 



P.B.I to 



-l)*'5^B , 



P.B.I to 



\Pa,Pb} 



(A. 17) 



P.B.I to 



0. 



Also for two generic functions, F(p, q) and G{p, q), which do not have explicit time 
dependence, we get up to the power, k. 



^([^,GWIto) 



G 



P.B. 



+ 



to 



P.B. 



to 



dF 
d^' 



G 



+ 



P.B. 



J 



P.B. 



(A.18) 



[F,Ht}p,j,,,g] + \f,[G,Ht}pA 

) P.B. L } 



P.B. 



This shows that the relations, (A. 17), actually hold up to the power, /c + 1, completing 
the proof. □ 



• Eq.(5.10) 
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The first equality in (5.10) follows from the algebraic identity (2.34) on the tangent space 
variables (g, q, t) . Firstly, one considers (5.7) for A — rh, 

d{5fidL_ ^ dm 

Qqm QqB Qqm ' ^ ■ ^ 

Secondly, one takes the partial derivative of each side of (A. 19) with respect to , 

QqA Qqm QqB y ' Qqm QqA QqB QqA Qqm ' ^ ■ ^ 

Thirdly, the partial derivative of each side of (5.6) reads explicitly as 

dq'^ dq"' dq^ dq"' dq^ ' 
where we made use of (5.5). Fourthly, making use of (A.20) in (A.21) leads to (5.10). 

The second equality in (5.10) holds because of the integrability condition (2.37). Proceed- 
ing step by step, one may start by using the chain rule in order to show the identity 

. ^^!^ ^ . l)(#.+#a)#.^ ( + d±d{5f)\ d-L 

y > Qqm dq^dq^dq^ ^ ' dq^ y dq"^ dq"^ dq^ ) dq^df 

(A.22) 

Then, the relation (2.37) is used to exchange some indices in Eq.(A.22) as follows 

^ ' y dq^ dq^ dq^ ) dq^dq^ 

^ ^ dq- \dq^dq^^ dq^dq^dq^ I ' ^' ^ 

Finally, one observes that the sum of terms in the parenthesis of Eq.(A.23) vanishes since 

O 'rh r\ f-\ • ID O 'rii ' \ * / 



dqriiQqB Qqm QqbQqB dq 
due to (5.5). The set of Eqs.(A.22)-(A.24) implies that 

'd'^L{q,q,t)\ _ dpeiq, Pa, Q"^, t) 

which ends the proof of Eq.(5.10) due to (3.7). 
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B Grassmann algebra 



In principle, in order to be able to discuss rigorously a generic dynamical system (i.e. which 
contains both bosons and fermions), one needs to introduce the "Grassmann algebra" Af^ which 
is generated by the N anti-commuting Grassmann variables [10], 



They generate the following basis for the Grassmann algebra 

1 } "body" 



(B.l) 



> "soul" 



(B.2) 



CC^ « < /3 



which has the dimension, dim(A^) = 2 while N can be infinity . 

Any quantity in the Grassmann algebra, Aj^, can be expressed as an expansion in terms of the 
above basis over the field IR of real numbers field (or the field C of complex numbers). 



N 



1 



X 



nl 



(B.3) 



n=l 



where xq and Xaia2-an real (or complex numbers) carrying totally anti-symmetric indices. 
Naturally, the bosons allow the expansion of even n's only, while fermions allow only odd n's. 

It is crucial to note that if and only ifxo 7^ 0, the inverse, x~^, exists. 

It is convenient to rename the elements in the basis with a given ordering as 

{Z^ 1< J<2^} = {l, r, C'C'---C^} , (B.4) 
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and to write 

2JV 



X 

j=i 

We also introduce the following notation to pick up the real or complex number coefficient. 



[x]j = xj. (B.6) 



In terms of the Grassmann algebra, the Lagrangian is a bosonic variable, but not necessarily a 
pure 'body'. Furthermore the actual dynamical variables are the real numbers, [g^] j, leading to 
a much bigger phase space. 

In deriving the equations of motion for a Lagrangian, what we actually encounter is the ex- 
pression, 

J dtSC{qn,t) = J dtS[q%Z''-^C{q^,t). (B.7) 

This implies that for the bosons, by considering especially the variations of the pure body, the 
equations of motion can be indeed collectively expressed as usual (2.7), but the equations of 
motion for the fermions should be refined to hold in a weaker form, 

= , J 7^ 2^ , (B.8) 

J 

when N is odd. Namely for the fermions, the usual equation of motion is true except the highest 
order in 'soul'. However this subtle issue can be neglected either by imposing the missing equa- 
tion for J = 2^ by hand, or by letting iV — > oo. 

The purpose of the present subsection was to provide a rigorous way to analyze the dynam- 
ical systems containing both bosons and fermions. Nevertheless, in practice we will favor the 
Lagrangian systems which do not require any explicit use of the basis for the Grassmann algebra, 
especially when they are transformed into Hamiltonian form for constrained systems. 



Mm, t) 
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C Basics on supermatrices 

A generic (ni + 712) x (mi + 777,2) supermatrix, M, over a Grassmann algebra, Aj^, 
(see Section B), is of the form, 

M^( ^"i^'"^ \ ^ (Qi) 

where A, B are bosonic and ^, © are fermionic. 



The complex conjugation, transpose, and the Hermitian conjugation read respectively [10], 

M* = ( Q I ( "J ) (M*U = (-l)#''(#« + #^)(M.,)*, 



(C.2) 



(C.3) 



Note that 

(M*)* = M, (Mt)t = M, (M^)t = M*, 

(M1M2)* = M*M* , {M1M2Y = MjMl , (iV/iM2)t = M^m/ . 

However, 

{M^f ^ M , (M*)'^ ^ {M^y , (Mt)^ 7^ M* , ete. (C.4) 



In particular, a real supermatrix is of the generic form, 

M = M* = Q ) ' = , (C.5) 

where every variable is real, A ^ A*, B ^ B*,"^ ^ 6 = O*. 
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For the (rii + 712) x {rii + ^2) square supermatrix, M, 



M= ( ^'^i^'^i ^"^^"^ ) , (C.6) 



the inverse can be expressed as 



^ \ -B-'Q{A-^B-'Q)-' {B-QA-'^)-' J' ^^''^ 
where we may write 

00 

{A - ^B-^Q)-^ = A-'^ + J2 {A-^'^B-^eyA-^ . (C.8) 

Note that due to the fermionic property of ^, 6, the power series terminates at p < nin2 + 1. 

The supertrace and the superdeterminant of M are defined as [10]^° 

strM = tr^-trB, (C.9) 

sdetM = det(^ - ^B-^Q)/detB = detA/det{B - BA"^^) . (C.IO) 

From Eq.(C.7), 

sdet M ^0 det AdetB 
is the necessary and sufficient condition for the existence of M~^. 

The supertrace and the superdeterminant have the properties, 

str (M1M2) = str (M2M1) , sdet (M1M2) = sdet Mi sdet M2 . (C. 1 1) 



For a generic n x n bosonic matrix, A, in a similar fashion to above, decomposing it into the 
'body' and 'soul' (see Section B), 

^ = Aody+^oul, (C.12) 



^"The last equality comes from det(l - A'^'i/B'^Q) = det~i(l - B'^OA'^*), which can be shown 
using det(l - X) = exp (- X)^! ^trX^) , and observing tr (^-1*5-19)?' = -tr {B-^QA-^'^)p . 
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we have 

^ ^ ^ ^body + X/ (~^body"^s°"l)''"^body • (C.13) 
p=l 



TTjms, ^4 ^ if and only «/^body 
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D Lemmas on the canonical transformations of supermatrices 



In this appendix, we do not explicitly state which entries in the supermatrices which are Grass- 
mann even or odd. The way the supermatrices have been written is supposed to be self-explanatory. 

Fact 1. 

For any n x m bosonic matrix M over a Grassmann algebra Aj^ (see Section B), 

M ^ (viV2 ■ ■ ■ t'mj , V.j = {Vij V2j ■ ■ ■ Vnj)'^ , (D.l) 

there exists an m x m nondegenerate matrix Q satisfying 

MQ = (wi W2 ■ ■ ■ Wk Si S2 ■ ■■ Sm-k^ , (D.2) 

where all vector (i.e. columns) Wi are orthogonal to each other and each body of them is nonzero, 
while the s^'s are pure souls or zero. 

Proof: To show this, one needs to separate the vectors Vi into two groups: pure soul ones 
and other ones with nontrivial bodies. Then one only needs to orthogonalize^^ the latter. □ 



Fact 2. 

For any n x m bosonic matrix M over a Grassmann algebra A^, there exists two nondegenerate 
matrices, P and P', which transform M into the canonical form, 

P'MP = ( ^'^^'^ Ofex(m-fc) \ ^ 

\ ^{n—k)xk S(n— fc)x(m— fe) / 

where Afexfc isak x k nondegenerate diagonal matrix, 

Afexfe = diag(Ai,A2,---,Afe), (Ai)body 7^ , l<i<k, (D.4) 

of rank k so that k < m and k < n, and the (n — k) x {m — k) matrix S(n-k)x(m-k) is a pure 
soul. 

^^Notethat, since the body is nonzero, the inverse of the scalar product («;tu;j)~^ exists, and the orthogonaUzation 
can be done by a generalization of the Gram-Schmidt procedure for the graded case. 
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Proof: From Fact 1, we construct P' out of the orthogonal vectors {ty|, • • • , w\} and 
their complementary vectors to get 

P'MP =i I f V (D.5) 



Now we only need to take one more step for the final result, 

A s' \( I -A-h' \ _ ( ^ 

□ 



(D.6) 



Fact 3. - Corollary 

For any nondegenerate n x n bosonic matrix M over a Grassmann algebra A^, meaning 

detMbodyv^O, (D.7) 
there exist two nondegenerate matrices P and P' which transforms M into the identity, 

P'MP = 1 . (D.8) 



Proof: The proof is straightforward from Fact 2 and its Proof, since the matrix is non- 
degenerate and n = m = k.n 



Fact 4. 

For any nondegenerate (ni + 712) x (711 + 712) square supermatrix M over a Grassmann algebra, 

M=(^'^ detAdet^T^O, (D.9) 

there exists two nondegenerate supermatrices P and P' satisfying 

P'MP^l. (D.IO) 
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Proof: After the transformation. 



(-x^^-^ b){o s-xV^)' 

we only need to apply Fact 3. □ 



Fact 5. 

For any nondegenerate, bosonic, real, symmetric or anti-symmetric matrix A± over a Grassmann 
algebra A^, 

A± = Ao + A,oni , = ^± , 

(D.12) 

(det A±)body = (det ^o)body 7^ , Al^ ±A± , 

there exists a nondegenerate real matrix P, the transformation induced by which, removes the 
pure soul Agoui completely, 

PA^P'^^Ao. (D.13) 

Remark: The point of Fact 5 is the removal or addition of any soul to the original matrix, A±, 
via the insertion between two appropriately chosen real matrices (D.13). 

Proof: We present explicitly the real transformation, 

oo 

P = P* = 1 + ^ a„(Ao,iAo-i)" , (D.14) 

n=l 

where the coefficients are given by a recurrence relation [15], with ao — 0, ai — 

n 

On+i = -On - I X] ajidn+i-j + ^n-j) , for n > 1 . (D.15) 
Due to the Grassmannian property, the sum in (D. 14) terminates at a finite order. □ 

Fact 6. 

For any bosonic, real, symmetric or anti- symmetric matrix, A±, over a Grassmann algebra, A^, 

A±^Al, Al^±A±, (D.16) 
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there exists a nondegenerate real matrix, P, which transforms A into the canonical form. 



PA^P^^ (^0^ °^ J , P^P\ (D.17) 

where the bosonic matrix, h± — ±6^, is nondegenerate, (det 6±)body 7^ 0, while s — is a 
pure soul or zero. 

Proof: First, using a real orthonormal matrix O, one can transform the 'body' of A into the 

h 



canonical form ( ^ g ) ' which gives 



T _ [ b + Si S2 



OA^O' = ( - j , (D.18) 

where h = ±6^ is nondegenerate and real, while si = ±sf , S2, S3 = isg are all real and 
pure souls. We further transform it as 

1 \ / 6 + S2 W 1 -{b + si)-h2 

To complete the proof, we only need to apply Fact 5 to 6 + si. □ 



Fact 7. 

For a generic (ni + 712) x (mi + 777,2) supermatrix M over a Grassmann algebra Aj^, 



M 



A \I/ 

®n2Xmi -^712 xm2 



(D.20) 



there exists two nondegenerate supermatrices P and P' which transforms it into the canonical 
form. 



P'MP^ 
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1 











\ 







Si 

















1 







\ 





X 





S2 


/ 



(D.21) 
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where Si, S2 are bosonic pure soul, and ijj, x are fermionic. The partition of the canonical form 
reads. 



[ki + (ni - ki) + k2 + (^2 - k2)] x [ki + (mi - ki) + k2 + (ma - /cs)] , (D.22) 

where ki, k2 are respectively the ranks of the bosonic matrices A, B. 

into the canonical form, in order to 



Proof: From Fact 2, we can first transform the xmi 
put M into the form 



1 





^1 





Si 




Xi 


X2 


B 



and to further have 




1 





^1 \ 







-^1 





Si 


^2 
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1 





Xi 


X2 









1 



Now we apply Fact 3 to (S — XiV'i) to get 



/ 1 



X3 



\ 



1 




V X4 

Finally, one completes the proof by the equality. 





S2 J 



1 1 











1 
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V 
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1 
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1 / 



/ 1 

Si 

X3 

V X4 





1 
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S2 J 



( 1 




V 






V 



\ 





X4 







S2 7 




(D.23) 





^2 

-B - XiV'i 
(D.24) 



(D.25) 



\ 

1 

-X3 1 

1/ 



(D.26) 



□ 
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Remark: Note also that Fact 3 follows as a corollary too. 



Fact 8. 

Consider a (n_ + n+) x (n_ + n+) anti-Hermitian supermatrix, over a Grassmann algebra A^, 
which has the symmetry property, Q,ab — —{—^)*"'*''^ba, or equivalently 



It is of the general form. 



n 



-1 



(D.27) 



(D.28) 



where every variable is real, A± = A^, ^ — and A± = ±A^. 

We note that the anti-Hermiticity and symmetry properties (D.27) are preserved under the 
transformations by a real supermatrix, (C.5), 



smce 



(D.29) 
(D.30) 



The claim is that there exists a nondegenerate real supermatrix, L — L*, which transforms 
into the following canonical form. 



(D.31) 





( 











\ 







s_ 































is+ 


/ 


where all the variables are real, 6± = 






i, = 


-- ip* 





matrices, (det 6±)body 7^ ; s-t are pure souls ; and b-i- = ±6^, s± = . 
The partition reads 

r 1 ^ 

A;_ + (n_ - A;_) + k+ + (n+ - k+) , 
where kzf are respectively the ranks of the bosonic matrices, A^, so that (&±)a;±xA;±- 



(D.32) 
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Proof: 



From Fact 6, we transform A_ into the canonical form, in such a way that 







Vi 



+ 



where = 



P 
1 



. We take it further to 




6_ 





















1 








: 








1 





^2 



Now apply Fac? 6 to + ip^hj-ipi) to get 






V 





s'_ 

T 

-X 




X 



•+ 



\ 




Finally, to complete the proof, we only need to take the following transformation. 



/I 





M 




( h. 





1 




' 
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/ 1 
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-^b+'x'' 
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1/ 



□ 
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E A paradigmatic example 

As an illustration of the general case, discussed in the core of the text, a most simple case is 
presented here: a Lagrangian 

L(g^g^) = ^(g^)-V-(g^), (E.l) 

which 

1. is a function on a bosonic tangent space of finite dimension A^, 

2. is of homogeneity degree equal to two (i.e. the system is free), 

3. leads to a positive definite energy, and 

4. does not lead to tertiary constraints. 

Hopefully, this example combines three virtues: (i) its simplicity should allow to displace the 

focus from the technical onto the conceptual, (ii) it includes both cases of first and second class 
constraints, and (iii) it provides the starting point of usual perturbative expansion, so it is not 
merely academical. 

The Lagrangian is assumed to be quadratic, therefore the kinetic energy T{q) — ^TAsq'^Q^ 
and the potential energy V{q) = \VabQ^Q^ are both quadratic forms. The Lagrangian (E.l) 
leads to a conserved energy equal to 

E{q,q)^T{q) + V{q), (E.2) 

which is positive definite, E{q, q) > 0, if and only if the kinetic and potential energy are sepa- 
rately positive definite: T(q) > ,V{q) > . Without loss of generality, one may assume that the 
variables q"^ have been 'rotated' so that the symmetric matrix Tab — d'^L/dq^dq^ is diagonal: 

a 

where the index a corresponds to the — strictly positive eigenvalues Ma > , while the 

index m corresponds to the remaining M. zero eigenvalues. The corresponding momentas are 
respectively given by Pa = MaCf (no sum on the index a!) andp^ = .^^ Therefore, one gets A4. 
primary constraints </)m(g, p, t) = and the primary constraint surface V is the hyperplanep^ = 

^•^In this example, notice that the distinction between hatted and unhatted indices is not necessary. 
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of codimension A4 embedded in the 2j\/'-dimensional phase space. The canonical Hamiltonian 
(3.24) reads 

H(q^p„t)\^ = _ ij] W?-. (E.4) 

a " A,B 

The total Hamiltonian (3.30) is thus given by 

HTiq^,PB, u^, t)^ H + u^pm. (E.5) 

The time evolution of through the Poisson bracket with the total Hamiltonian leads to the 
equality = {g™, Ht}p.b. = u™- . The one-to-one maps (3.34) can be seen explicitly in the 
present case since Pa — MaQ"' (no summation) and — q"^ . 

The preservation (4.18) of the primary constraints under the time evolution leads to the sec- 
ondary constraints 

t) = Vmraq"^ + \Y. ^rnnq'' + \ ^rna^" • (E.6) 

Now the point is that the preservation (4.43) of the secondary constraints under the time evolution 
would leads to new, i.e. tertiary constraints if some entries Vma were non- vanishing. Therefore, in 
the particular example we are considering, one assumes that the potential does not include mixed 
terms: 

V{q\ (r)-\Y. + ^ E ^-n?"^g" ■ (E.7) 

For that reason, one may perform a rotation in the plane of the variables g"/ y/Ma in order to 
make the symmetric matrix Vat diagonal without modifying the kinetic energy. Without loss of 
generality, the symmetric matrix Vmn also be assumed to be diagonal. Therefore, 

V(q% q\ q^, q^) = ^ ^(a;^)2(g^2 ^ 1 ^(^-)2 . (g.g) 

a m 

where the 'barred' indices correspond to the strictly positive eigenvalues while the 'Greek' indices 
correspond to the vanishing eigenvalues of the matrices Vgh and Vmn ■ The secondary constraints 
(E.6) become simply 0^(g, p, t) = q™- . There are no tertiary constraints because the preservation 
(4.44) of the secondary constraints under the time evolution only leads to the fact that the La- 
grange multipliers — Q , while the it'* 's can be arbitrary functions of time, which signals the 
presence of some gauge freedom. The constraint surface V is thus the hyperplane defined by the 
system = p^: = g" = . It is straightforward to check that the constraints are (primary) first 
class constraints and that the set {pm, q^} contains all the second class constraints. Notice that, in 
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this example, the total and extended Hamiltonians are identical since there are no secondary first 

class constraints. 

Under the sole hypotheses stated above, the Lagrangian and the Hamiltonian have been de- 
composed into a sum of four pieces^^ 

^ _ i^^QQ _|_ ^harmonic _|_ ^^Ist _|_ ^2nd 

^ _ ^free _|_ ^harmonic _|_ ^ 1st _|_ ^2nd (E 9) 

where each piece corresponds to one of the following four distinct physical cases: 

• Free particles: 

a 

o-free/^ x _ \^ 

a 

• Harmonic oscillators: 

^harmonic(^. = ^ [m,(^^)^ - ■ 

a 

^harmonic(^.^) = ^T^i^ + M'i^^] ■ 

a 



• First class variables: 

^Ist/^ „,u\ \ ^ „.tt„ ijlst 



List = 0, 



Second class variables: 



m 



2 

m 



^^The present "paradigmatic" example has been inspired from the two examples given in the section 1.6.2 of [7] 
which correspond to L^^*^ and L^^'^. The straightforward quantization procedure for these two cases is respectively 
carried on in the sections 13.1.1 and 13.1.2. 
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